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ABSTRACT 


A semi-analytical solution to the problem of the motion 
of a satellite of the moon is presented. The theory is de- 
veloped to third order, where first order is 10'^. Pertur- 
bative effects which are considered include those due to the 
attraction of the moon, earth, and sun, the non- sphericity of 
the moon's gravitational field, the oblateness of the earth, 
coupling of lower-order terms, solar radiation pressure, 
and physical libration. Short-period terms and those with 
the period of the moon's longitude are produced by means 
of von Zeipel's method; it is proposed to obtain the secular 
perturbations, and those depending only on the argument of 
perilune, by numerical integration of the equations of 
motion. 
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A SEMI-ANALYTIC THEORY FOR THE MOTION OF 
A LUNAR SATELLITE 


by 

Giorgio E. O. Giacaglia* 

Yale University 
and 

James P. Murphy and Theodore L. Felsentreger 
Goddard Space Flight Center 


INTRODUCTION 

The motion of a satellite of the moon, or lunar orbiter, is analyzed. The solution is developed 
in powers of lO'^-i.e., first order is 10 ^ second order is 10'% etc. The Hamiltonian for the 
"main problem" consists of zero-, first-, and second-order quantities. The higher-order Ham- 
iltonian is of third order and consists of two parts, the first containing terms generated by coupling 
of lower-order terms and the second consisting of terms added by considering further perturbing 
forces, such as solar radiation pressure, physical libration, non-sphericity of the earth’s potential 
field, the attraction of the sun, etc. Additional terms are produced by considering also the ec- 
centricity and inclination of the moon's orbit. 

In order to retain the relative orders of the disturbing forces, it is necessary to restrict the 
semi-major axis of the orbiter to about four moon radii or less. Furthermore, certain restrictions 
must be made on the eccentricity and inclination in order not to invalidate the solution. There- 
fore, the following assumptions on the semi-major axis, eccentricity, and inclination are made: 


a ^ 4 moon radii 
.01 < e < .75 
sin I > . 01 . 

The small parameter of first order is n*, which is the mean motion of the moon's mean longi- 
tude. Recent determinations of the spherical harmonics of the moon (see References 1, 2, and 3) 
indicate that the second degree zonal and sectorial harmonic coefficients, together with the third. 


*Now at the University of Sao Paulo, Brazil. 
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fourth, and fifth degree zonal harmonic coefficients, are all of about the same order of magnitude- 
namely,10"‘. Therefore, the small parameters of second order are Jj, J 22 , J3J J4> Js» ^.nd (ri(j/n)2, 
where J2 and J 22 define the principal part of the oblateness of the moon and /n)^ is the square 
of the ratio of the mean motion of the moon to that of the orbiter. The quantities J3, J4, and J 5 are 
higher degree zonal harmonic parameters of the moon. The small parameters of third order are 
j,(ng/n)S (n^/n)2, (n^/n)= sin (i^/2), K/n)^ e^, a, {n^/n)% and anjn. Here, the parameter J2 
is the principal part of the oblateness of the earth. The quantity (n^/n)^ is the square of the ratio 
of the mean motion of the earth to that of the orbiter. The fact that the moon s orbital plane is in- 
clined to its equator and the fact that the moon's orbit about the earth is elliptical give rise to the 
two small parameters of third order sm(ig/2) and respectively. The radiation 

pressure gives rise to the third order parameter cr, and (n^/n)^ is the cube of the ratio of the mean 
motion of the moon to that of the orbiter. Finally, is the correction due to physical libration. 

The longest meridian of the moon contains the line joining the centers of mass of the earth 
and the moon. The right-handed, rotating, selenocentric coordinate systems adopted for this prob- 
lem will then be as follows: The z-axis is the rotational axis of the moon, and the xy-plane is the 
moon's equatorial plane. The x-axis passes through the moon's longest meridian, and is assumed 
to rotate with the motion n* . 

For a semi-major axis of 4 moon radii, 

i^)‘ - 

(v) - 2.7 « 10 -‘ 

so that the oblateness terms in the lunar potential and the earth perturbations are both of about the 
same order. The largest oblateness term of the earth is then of the order of 10"® . The perturba- 
tion of the sun is of third order. 

Another perturbation to be considered is the effective radiation pressure of the sun, whose 
strength is about 1 x 10 ^ dyne/cm^. If the area— mass ratio of the orbiter is 1.5 x 10 cm /g, 
then the disturbing acceleration due to radiation pressure is also of third order. 


EQUATIONS OF MOTION 

The first step is to determine the equations of motion for the gravitational fields of the moon 
as a primary and the earth and sun as perturbations. 

The following notation will be used: 

subscript 0 : moon 
subscript 1 : orbiter 
subscript 2 : earth 
subscript 3 : sun. 
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In an inertial system, the equations of motion are 


m. p. - grad^ U (j - 0, 1, 2, 3) , 


where 


u = k" 




( 1 ) 


and p. is the radius vector of any one of the four bodies. If 77., are the rectangular inertial 
coordinates of one of the bodies, then the equations of motion can also be written 


vu 

wr a ^ O' 2 . 3) . (2) 

Similar expressions hold for 77. and ? . . 

It is now convenient to refer the orbiter to a moon-centered system, the moon to an earth- 
centered system, and the sun to a system whose origin is at the center of mass of the earth-moon 
system, designated by the subscript G. Therefore, 

- -fo ' 


and 


X 


3 


^3 - 


■"o ^o'""’2^2 


( 3 ) 


Similar expressions hold for and z.. 

The equations of motion must be transformed accordingly. The partials in Equation 2 are then 
computed with respect to the new variables by making use of 


(9U 



au 

(k = 0. 1. 2. 3 ) . 

j 


j = 0 
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Ill 


I I I II I mill ■ 


It follows that 


6V 

dU 


-iLL 

(9^2 


au 

3x„ 


au 


au _ "'o au 

axj “ axj 


au 

ax„ 


"2 au 


"0 ■^'"2 ^^3 


and 


Since 


au 

a^3 


au 

<9xo 


- ■ ^0 


(4) 


it follows that 


.. _ /j_ J_\ ^ _ J_ ^ -^U . 

^1 " ^nij m^J axj mg ax^ +mj axj '■ 

The force function U must now be expressed in terms of the new coordinates. We have 


o\ = (^0-^2)' + ••• " ’‘o' ^ ••• " 

o'a = (^3 -^o)' " ••• ^ 


,2+2+2 
‘1 ^1 1 


x/ + y/ + z/ + 


^3 nig + mj ’‘o 




2m, 


) - + 1 ^; (’‘o ’‘3 ■ Voya ^ ^0 ^3) 


2 m„ 


^ 3 " ^ lo ' o '^'" 2 y "’ o '^"’2 '^0 ’ "^3 • 




1 * 0 0 
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and 



(^2 ^3) ^ + • • • 


‘3 nijj +m 2 ^ 




r,'' + 


mo +mj 


2m„ 


r " + 


"’0 "'2 


Since U contains the inverses of these radii, note that 


(6) 



The angles used are shown in Figure 1 . The angle S03 may be expressed in terms of So'3, as 
follows: Since 


and 


’”0 "'2 
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Figure 1— Relative positions of the moon, 
orbiter, earth, and sun. 


Finally, 




= 77 - 


f3' ™0 + "'2 


Certainly, to third order, - r^, so that 


COS 


cos Sq'j - 



(8) 


On the other hand, may be replaced by 180°— where S/p is the selenocentric elongation 
between the earth and the orbiter. 


OBLATENESS TERMS 

If the plane of reference is the lunar equatorial plane, then the disturbing force per unit mass 
may be written as 


7 2— i2_ivi/ 


j)cosm(\-A„J . 


The only terms to be included are 522} J20 " “J2J J30 ^ "Js’ J40 " " J5 • Therefore, 


u 


OBL . 





J P (sin f}) 

n n ' ^ 



J22 P22 cos 2(\ 


(9) 


The angle /3 is the latitude of the orbiter with respect to the equator of the moon, X' its longitude 
reckoned from any fixed direction, and X22 the longitude of the moon's longest meridian from the 
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same fixed direction. Note that and can be expressed in terms of the coordinates x, , yj, z 
of the orbiter. However, will contain the time explicitly, since 


^22 ■ ^22 (0) + t , 


where is the frequency of rotation of the moon around its axis. Further, if we neglect 
the physical libration of the moon, the longest meridian is always pointing toward the earth 
and - rig. 


CANONICAL EQUATIONS ANO GRAVITATIONAL TERMS 


Let us choose as canonical variables the Delaunay set 


L 

G 

H 


/mq 3 . I = mean anomaly, 

L n - , a = argument of pericenter, 

G cos I , i; = longitude of ascending node. 


where a, e, and I are the semi-major axis, eccentricity, and inclination, respectively, and where 
/'o " "ig, in which k is the Gaussian constant and is the mass of the moon. Then the equations 
of motion become 


L 


'21 
ri I 
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H 


2Z 
ao - 


<9F f?F 

-9L ■ " ' “ r(G • 


n 


2£ 

■<9H ' 


where 


•0 

2L2 


+ U.. 


In Equation 11, is a function which has to satisfy the condition 

•• - 

’‘l ^^GRAV. ^ ‘ 


(10) 


( 11 ) 


( 12 ) 
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The force fxmctioii u depends on the new variables x through Equations 6. Then, 


3U 

5x„ 


= k2 


-nig 3 nig 

02 


"o +">2 ’'o / ">2 \ ^1 ^0 


^1 ^3 itIq + nij 


. 3 
1 3 


V^O ■"'"2/ 


irlg + m 2 


d 

dx. 




. 3 

02 


+ k = 


'”0 ™2 '"3 

^3 • ^1 

, mom/ m 3 

^0 

• ^1 

+ k2 

1 


mg +mj 


^ / \ 9 

(mo 

r 

3 

03 

"’T 2 ri2 


m 

. 1 7 ^ 

m 2 m 3 1 

^ mg mj m 3 

^3 

* 

1 2 

mg^m^mo 

^0 • '^1 

+ k ^ 

mg 

+ ™2 •'13 

~ mg + 

r 

3 

23 

~ k 

(mg +m^}^ 

^2% 


(13) 


and 


3U 

i9x. 


k" I -nigmj 


^3 mg + ’‘o 


^^3 m„ + m„ ^0 


0 "2 


’'s mg + m^ ^0 


<9x, ■ "’o'"3 , 3 


3 ‘ 1 , ">0"'2'"3 '•0 • "-j 

<■ k^ 


m + m 3 

” ^ *^03 


’^"'"’l"’3 777 - ’^^'"2™3 ^3 


3 ’ ^ 1 "^0 "’2 ”’3 0 ’ 1 


m +m 3 

u z 1^3 


(14) 


Then, substitution of Equations 13 and 14 into Equation 5 yields 


”1 ^ ™0 I? 


[tIq rrij 2 2 

FT + k 2 m,m 3 — 


mg 5Xj 


’^0 *^1 ,"' o "’ 2''’3 ’^ 3 ’^! 

,2 „ „ + k2 


1- m 


0 2 ^ 3 


"'o + "’2 rA 


- k^ 


^ ■■■ 7'F + -"T2 F77 


(mg +m 2 )- 


03 
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j 2 

m m 2 ni^ 

1 

m ^ m A m o r « * r 

,9 'J >5 -i 

1 . , "’o^"’2">3 

1^0 • 

k 

mg + m 2 ! 

^13 

mg+m2 J .233 

k . « 2 

(™o +"> 2 ) 

^2^3 ) 


1 ^ 

, 0 

1^3 ■ 1^1 

mg m 2 m^ r p • r j 


ITlo 

+ m2 i5xj 

- 

"'o'"3 ,3 + 

^03 

("lo + " 12 ) Fg33 





_ — ♦ 

_ „ _ - . -n 



- k^mjiTij — “ ^2 11,2013 -yj 


Fj • Tj trig tllj ITI 3 Tg • Fj 


5x, 


"2 -3 

r“ " 7 : 


'‘3*1 ‘3 _ 3 

03 , 


Therefore, upon comparing Equations 12 and 15 we find 


mj oij O’ 3 


''o • ■'i 

- - + k-'mj ^ 


m. 


where 


F, - F 


3 mg + nij I'd 


and 


series in Legendre polynomials, 


F = series in Legendre polynomials. 


The term k^mj/Fg^ can, of course, be neglected. Finally, from Equations 9, 11, and 16, 


F = 


Mo M 2 M 3 

2L^ OJ2 Oj3 


■ 0 • I'l M3 Ti • F3 


-£(^) (7 

n= 2 


7y J 22 P 22 (sin/3) cos 2(\' 


• (15) 


(16) 


• (17) 
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^r> 


Now 


- rg Fj cos Sj’o 


1 S/o 


and 


Then, 


<» . « 

1 _ 1 '‘i ,1 / , \ 

r — “ T“ ^ — o + — / P (cosS.fi) I — 

>^12 "^0 r 2 10 r„ p V 10/ \^F(, 


2 ^ 0 '■ I 

~~ ^ ^2 

1 o 


/^2 Mj /rjN^*' 

^ r„ ^ r,2^ irj ?p 


P=2 


and since independent of the position of the orbiter, Equation 17 becomes 


1-^0 ^2 \ ' j " 1 \ , ^3 ‘ I ‘ 3 

2L^ * I'd Z_, Uo) (^°^^io) + Fj 3 ■ 1^3 rA’ 


p=2 


5 /r> \ n 

y~^ (^j j 


^\T;I JnPn("i"^')^ i^T7j J22P22^^i'^'-^)^°^2(\'-\, 

n - 2 

A further simplification can be made if one considers the following expansion: 


Then, 




^1 ^ Pq3 2Ti ^03 ‘^°®^13 


^0^3 


1 + 1 ) - 2 cosSj" 

^n-s 


Pl3 ^ 103 ^(^ 03 ) 


1 Pi^osS/'j 1 ^ /rjNP 


’’03 2Z ('' 03 ) ^P ■ 


p=2 


) • ( 18 ) 
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Also. 


1 *^03 ^13 


Ti COS S/j 


. 3 
03 


Since is independent of x^, y^, z^, Equation 18 becomes 


^ (rj Pp (,^‘J (cosS;3) 

p=2 p=2 


r , 


2 ^( 77 ) JpPn^^in/^) ^ J3,P22(^i"/?)'^°^2(\'-\33) 

n= 2 


(19) 


From now on, the subscript 1 will be omitted. 

Including terms of second order, the Hamiltonian becomes 




2L^ 


77 ( 77 ) [P 2 * 77 P 3 j Pj (cosS/'j) 


JnPpC^i'^'^) * (-^) J22P32(sino)cos2(x'-\37 


( 20 ) 


THE ANGLES S AND S ^3 

In view of the form assigned to the oblateness terms, the plane of reference is the equator of 
the moon. Therefore, the next step is to express the angles Sj'g and Sj^ in terms of the orbital 
elements of the orbiter, the moon, the earth, and the sun, with respect to that plane. The geometry 
is shown in Figure 2. 

The explicit form of s/p and requires the solution of two spherical quadrangles. This is 
now done. 

If X, y, z and x^, y^, Zg, are the rectangular coordinates of the orbiter and the earth, respec- 
tively, and r and their selenocentric distances, then 
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b. The orbit of the moon around the earth is circular and the motion uniform. The error is 
proportional to the moon's eccentricity, or about 0.055. 

c. The sun's perturbations are negligible. The relative error for a moderately high satellite 
is about 0.05. 

d. The mean longitude of the earth, is equal to 


With these approximations, the precision of the disturbing function is not higher than lO"* . It 
will be called the "disturbing function of the main problem," and it is given by 

^ ^ ^0 f/M" r T 

^ ^ ~ " P 2 ; [- J2P2(sin/3)+J2,P,, Csin/3) cos2(\' -\J] 


I'M' . iM' 


JjPjCsin^)^ ■ (23) 


where 


>‘■0 


and where the angle s/g is as shown in Figure 4. Since ic = 0 and can be replaced by 

it then follows that 


cos S, 


cos V cos 


(fi - - sin V sin (0 - \^) cos I 


(24) 


The equations of motion are given by the ca- 
nonical set 


L 


G 


H 




dF 



SI ’ '■ SL 



SF 

So. 

" ■ <9G 

SF 


sn ' 

^ “ SE 


(25) 


Since the variables n and Xg, appear only in the 
combination where x^, = n^* t + const., the 

degree of freedom is reduced by one by choosing 



Figure 4 — Simplified selenocentric celestial sphere. 
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as a new variable 


h = n - . 

The Hamiltonian must be modified accordingly, for 


n - n„ 


<9 \ 


(i.e., upon substitution of F + n^* H in place of F^. 

The Hamiltonian is stiU time -dependent through Since the longest meridian is always 
pointing toward the earth, it is possible to choose the rotating system whose x-axis passes through 
this meridian. The final form of the Hamiltonian for the main problem is therefore 




F ^ 2^, + n,, H + , 


PjCs) + [' J2P2^sin^)+J22Pj2(sin^)cos2\] 

- (^) JjPjCsin/O- (^) J4P4 (sin/^) - (-^) J5P5 (sin/o| ■ ( 26 ) 


where 




X' 


The equations of motion are (using g in place of ^): 


SF 

L = ^ • 



H 


(9h 


<9F . _ ^ ■ _ _ ZL 

‘ <9L ’ ^ ^ ~ SG ' " <5H ■ 


(27) 


In the discussion that foUows, the solution of this system will be given using von Zeipel's method. 


DEVELOPMENT OF THE DISTURBING FUNCTION 

The notation used will be as shown in Figure 4. 

From Equation 24, 

cosSj'o = s = cos (v +n-X^) - 2 sin^ sin V sin (Xg, -n) ; 
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thus 


?2 (S) 



“sin f cos f 


+ sinfcosf 


cos ( 2 g+ 2 h) + sin^ (g + h) 

sin ( 2g + 2h) + 4 sin"* ^ sin^ h [cos^ g - cos^ f cos 2g 
sin 2 g] + 4sin^ 2 sinh[^sinfcosfcos(2g + h) 

-cosgsin(g + h) + cos^f sin( 2 g+h) 



Further, 


sin /) sin I sin v 

cos cos ( >■ - h ) 
cos sin (> - h) 


sin I sin ( f + g) , 
cos V . 
cos I s i n V , 


and therefore 


cos cos A " cos V cos h ~ cos I sin v sin h , 


from which 

D. N 3 1 3, - 1 

(sin ■■•) ysin^,'-^ -^sin"^Isin''(ftg)--j- 

Pj 2 ( sin /j) cos 2X " 6 cos^ a’ cos ^ \ - 3 cos^ /i 

" 6(^^cos^f+X^sin^f + 2^xsinfcosf) - 3(l-sin^Isin^v), 

5 3 5 2 

P 3 (sin/i) = -J sin^ /i - ^ sin /3 = sin® I sin^ ( f + g) - y sin I sin ( f + g) , 

P 4 (sin/3) - “g“ sin“* /i - “y sin®/3 + y = g sin‘’lsin‘*(f + g) - y sin®Isin®(f + g) '•■g'’ 
Pj (sin/3) = y sin® /3 -y sin® ^ + y sin^ = y sin® I sin®( f + g) - y sin® I sin®( f + g) + y sin 


• ( 28 ) 


sin( f + g) 

( 29 ) 
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where 


5 - cos g cos h - cos I sin g sin h 


and 

X ~ ~ sin g cos h - cos I cos g sin h . 

Therefore, writing for in Equation 26, we have 

o n.^ r ^ r 

F = ^ n * H + ^ ; -icos^ f cos ( 2g + 2h) + sin^ (g + h) - sin f cos f sin ( 2g + 2h) 

2l 2 e z £ 1^ 

+ 4 sin'* ■§ sin^ h [cos^ g - cos ^ f cos 2g + sin f cos f sin 2g] + 4 sin^ sin h [sin f cos f cos ( 2g + h) 


- cosgsin(g + h) + cos^fsin(2g + h)J 


^22 / A. 

^]| ■ 2« "c r 2 

y Jj + J 22 [e cos^ f t sin^ f + 2Cxsin f cos f) - 3 (l - sin^ I 


J sin^ I sin^ ( f + g) 


J 3 


g , , 3 , , 

J 4 

y sin® I sin® ( f + g) - J sin I sin ( f g) 

r® 


15 2 

- -y sin^Isin^(f + g) + ^ 


J 5 


sin"* I sin'* ( f + g) 


63 35 

-y sin® I sin® (f + g)“ "y sin®Isin®(f + g) 
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+ "y sin I sin ( f + g) 


(30) 


= F„ + F + F, 


where 


Fq = (0th order) , 

Fj = ng H (1st order) , 
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and 


Fj - F - (Fj+f„) (2nd order) . 


(31) 


ELIMINATION OF SHORT-PERIOO PERTURBATIONS 


The terms in F which depend on I will result in short-period perturbations. According to von 
Zeipel’s method, the elimination of these terms corresponds to the solution of the system 

Fo'(L') -- F, (L'). 

dSj ^F„ 

F,'(H') = F^(H') ^ ’ 

^2 ' ^2 ' r}h r)H' ' SI rlL' ' 2 / j ' (32) 

where the generating function of the transformation 

(L. G, H, /, g, h) - (L', G'. H', , g', h') 


is 


S L' / < G' g * H' h t Sj ^ Sj t • • • , 


(33) 


and the new Hamiltonian 


F ' ' F ' ^ F ' 

^0 i ^2 


(34) 


should be independent of / ' . 

In order to accomplish this last requirement, a particular solution is given by 


f ’ f 

*^2 2 s 


and 


0S2 < 9 Fq 


i)l ^L' 2p - 

where F^^ and F^p are, respectively, the parts of F^ independent of l and dependent on 1. 


(35) 
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If Sj is determined from this system, then the relations between old and new variables are 


dS^ ^ 

I' = I + ’ g' = g + 


h' = h + 


5S2 


L' + 


dS^ 

TT 


G = G' + 


aS 2 

<lg 


H = H' + 


< 9 S 2 


( 36 ) 


Next, the "secular" part of F2 is determined. By definition, 


•'0 


F^s ^ ^ 1 F2dZ , 


(37) 


Using the relations (e = eccentric anomaly) 


' ( 1 “ e cos E) dE , 


r cos f a(cosE"e) 


r sin f a iT- sin E 


and 


r " a( 1 - e cos E) , 


' a 2 (l 


1 * e cos f 

7 [l - e^) 


(38) 


it follows that 


2L 


+ n^H' + 


1 2 e 


"c 1/ G' 


1 t 3 


t 3 1 


G' ■ 


cos 2h ' 


+ 15 1 


G' 


2 


G' 


H'- 


cos 2 (r' ’ h' ) t (1 - —^1 cos 2g' + 2 U ~ 77 ) cos 2(g' - h' ) 


H' 




1 ^ (L'Y 

^ L' ® \G7 


Re J 2 (1 


H' 


3 I " J 


2 

ffi -> 22 


/ H '2 


cos 2h ' 


18 




Also, since 


L' 


ai' 


0 , 


then 


L' - const. 


The generating function Sj and the short-period perturbations are obtained next. Since 


''h 

rVL' 



from Equation 35 it follows that 


S 


2 




(40) 


The integration is carried out using E or f as independent variables, again making use of Equa- 
tions 38. As a result. 


’2 " 8 


1 

J 2 

[- 2 (l - 3 cos^ I 

) 

8 n' 

a'3 (i -e'2)3/2 

1 

J 22 

[e (sin2 I') 3^2, 1 

+ ( 1 - 

8 n' 

a'3 (l-e'2)3/2 
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+ 


+ 


+ 


1 J 3 sinX r n \ 

96 n' a'" (l-e'2)s/2 (l - 5 cos I ) B 3 j + (sin I ) 

®3,2 J 



512 3.5 2 ) 7/2 [-2(3-30cos='I + 35 cos'* I ) 

B 4 .I 



2 sin^ I ' (1 “ 7 cos^ I ') B^ 2 “ (sin'' I ') 3 J 




1 J 5 ^ ' r 

30720 a -6 2 ) 9/2 b°(l 14^0^21 421cos"l 

') ^5 

sin^ I ' (1 ” 9 cos^ I 

')B3,2-3(sin^I')B3j 

1 ^ ' ^ f r 

384 |l2[- 2 (1 - 3cosM') (2 + 3e'2) H 6 sin2 I' 

(2 + 3e'^) 

cos 2h + 30 e ' ^ s i n ^ 

I ' cos 2g 


15e'2 (1 - cos I')2 cos C2g - 2h) + 15e'2 (1 + cos I')^ cos (2g + 2h)] ■(£'-/) 

+ 4(1 -3 cos 2 I') Bg j ^ 6 (sin^ I') j ^ 3 (l - cos I ' ) ^ 3 + 3 (l + cos I ' ) ^ ^ j ' (41) 


where 

Bj j - e ' sin f ' + ( f ' - O • 

B^ 2 = 3e'sin(f't2g)<3sin(2f'-t2g’)-te'sin(3f'<2g). 

Bj 2 1 "" 2( f ' - / ) cos 2h ' e' sin ( f ' - 2h) ■ e' sin ( f ' ’ 2h) . 

Bj 2 2 " 3e' sin ( f ' + 2g - 2h) + 3 sin ( 2f ' ’ 2g - 2h) + e ' sin (3f ' + 2g - 2h) , 

Bjj 3 = 3e' sin ( f ' + 2g + 2h) + 3 sin (2f ' + 2g + 2h) + e' sin ( 3f ' + 2g + 2h) , 

B 3 j = 12 e' (f ' - O sin g + 3 e '2 cos (f ' - g) - 6(2 + e'2) cos ( f ' + g) - 6e' cos (2f ' + g) - e'2 cos (3f ' + g) , 

B 3 2 = lSe'2 cos (f' + 3g) + 30e' cos (-2f' + 3g) * 10(2 + e'2) cos (3f' + 3g) 

+ 15e' cos(4f' + 3g) + 3e‘^cos(5f' ^3g') , 

1 "" 6 (2 + 3 e' ^) ( f ' - i) + 9e' (4 + e' sin f ' + 9e' ^ sin 2f ' + e' ^ sin 3f ' , 
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, 2 60e' 2 ( f ' 1) cos 2g + lOe'^ sin (f ' - 2g) + 30e' (4 + e'^) sin ( f ' + 2g) + 20 (2 + 3e ' ^) sin (2f ' + 2g) 

+ lOe' (4 + e'2) sin (3f ' + 2g) + lSe'2 sin (4f ' + 2g) + 2e'^ sin (5f ' + 2g) , 

^4,3 = 3Se'3 sin (f ' + 4g) + lOSe'^ sin (2f ' + 4g) + 35e' (4 + e ' ^) sin (3f ' + 4g) + 35 (2 + 3e ' ^) sin ( 4f ' + 4g) 

+ 21e' (4 + e'2} sin (5f ' + 4g) + 35e'^ sin (6f ' + 4g) + 5e'^ sin (7f ' + 4g) , 

Bj 1 = 120e' (4 + 3e'2) (f ' - Z) sin g - 30(8 + 24e'2 + 30"*) cos (f ' + g) + 60e'2 (6 + e'2) cos (f ' - g) 

-60e' (4+3e'^)cos(2f' + g) + 60e ' ^ cos (2f' - g) 

- 20e' 2 (6 + e'2) cos (3f ' + g) + 5e ' •* cos ( 3f ' - g) - 30e ' ^ cos ( 4f ' + g) - 3e cos ( 5f ' + g) , 

^5.2 ^ “ 840e'^ (f ' - / ) sin 3g + 420e' ^ (6 + e ' 2) cos ( f ' + 3g) - lOSe"* cos ( f ' - 3g) 

I 420e' (4 t 3e'2) cos (2f' + 3g) + 70 (8 < 24e ' ^ + 3e"') cos (3f ' ^ 3g) * 210e' (4 t 3e'^) cos (4f' + 3g) 

+ 84e'^ (6 e'2) cos (5f' + 3g) + 140e'^ cos (6f' + 3g) + 15e"* cos (7f ' + 3g) , 

, 3 ^ 315e"* cos ( f ' + Sg) + 1260e'^ cos (2f ' + 5g) + 420e'2 (6 + e ' 2) cos ( 3f ' + 5g) 

f 630e' (4 + 3e'2) cos (4f ' + 5g) + 126(8 + 24e'^ t 3e"') cos (Sf ' t 5g) + 420e' (4 + 3e'^) cos Cef' + 5g) 

+ 180e'2 (6 t e'2) cos (7f' ^ 5g) ♦ 315e'^ cos (8f' + 5g) + 350"* cos (9f' + 5g) , 

^ (4 + e'^) sin E' - 9e'2 sin 2E' + e'^ sin 3E' , 

^e,2 " - 9e' (4 + e' sin (E' - 2h) + 9e'2 sin (2E' - 2h) - e'^ sin (3E' - 2h) 

- 9e' (4 + e' 2) sin (E' + 2h) + 9e ' ^ sin ( 2E' + 2h ) - e ' ^ sin ( 3E' + 2h) 

- 15e' [(2 + e'2) - 2(l- sin (E' - 2g) + 3 [(2 + e ' 2) - 2 (l - e ' ^ ) 1^2 (l + e'^)] sin (2E' - 2g) 

-e' [(2-e'2) -2(l-e'2)'"2] sin ( 3E' - 2g) - 15e ' [(2 + e ' 2) + 2 (l - e ' = ) ' sin (E' + 2g) 

+ 3 1^(2 + e' 2) +2(1- e' 2)1/2 (l +e'2)j sin(2E' + 2g) - e ' [(2 - e ' 2) + 2 (l - e ' 2 )i /2j sin (3E' + 2g) , 
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Bg 3 = - 15e' 1^(2 + e'2) - 2(1 - gin _ 2g + 2h) + 3 [^(2 + e'2) 

- 2 (1 - e' 2 ) 1/2 (1 + e'2)] sin ( 2 E' - 2 g + 2 h) - e' 1^(2 - e' 2 ) - 2 (1 - e' 2 ) sin ( 3E' - 2g + 2h) 

- 15e' [^(2 + e' 2) + 2(1 - e' 2 ) 1 / 2 J sin (E* + 2g - 2h) + 3 j^(2 + e ' 2 ) + 2 (l - e'^) 122 (l + e ' 2)Jsin ( 2E' + 2g - 2h) 

- e' 1^(2 -e' 2 ) + 2 (1 - e' sin ( 3 E' + 2g - 2h) , 

and 

4 = - 15e' [^(2 + e'2) - 2(1 -e'2)i22j gin(E' - 2g- 2h) + 3 [(2 + e'2) 

- 2 (1 - e' 2 sin (2E' - 2g- 2h) - e' [^(2 - e'2) - 2 (1 - e' 2 )’/ 2 ] sin ( 3 E' - 2g - 2h) 

- 15e' [^(2 + e' + 2 (1 - e' sin (E' + 2g + 2h) + 3 j^(2 + e ' ^) + 2 (l - e ' ^ ^ ^ (l + e ' ^)j sin ( 2E' + 2g + 2h) 

-e' [( 2 -e' 2 ) + 2 (l-e' 2 )'^ 2 j sin^E- + 2 g + 2 h) . 

Thus, for the short-period terms, 


S = L' / + G' g + H' h + 83 

to second order. 

The short-period perturbations are given by 


/' = 


h' 


L = 


(9L' 

dC 

.98 

.9H' 


.98 
d I 


I + 


h + 


L' + 


dS^ 

dh' 

SG' 

.9S2 

.9H' 

dS^ 

yr 


c)S 

<?g 


G' + 


.982 


( 43 ) 
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and 


H 


as 

ah 


H' + 


as^ 


(44) 


It is easier to compute the partial derivatives with respect to the Keplerian elements a', e', i' and 
then compute them with respect to L', G' , H', as follows: 


asj 

^^2 da' ^^2 

de' 

L' <^^2 G '2 '^^2 

ah' 

da' ah' de' 


da' e' L'^ <9e' 


as^ 

as^ 

de' 

as^ 

ai' 

G' "^^2 H' 

<983 

dG' 

de' 

ac 

+ 

<91 ' 

dG' 

de* G'^sinl' 

-91 ' 


and 


<983 

-•98 3 

ai' 

1 

ah' 

-91 ' 

ah' 

G ' s i n I ' f9 1 ' 


(45) 


Furthermore, it is important to note that 


ar 

ae' 




sin f ' 


r9E' 

ae' 


s in E' , 


and 


<5n ' 3 n ' 

da' 2 a' (46) 


The partial derivative with respect to I does not need to be computed again, since 


as^ 

'TT 



T (F-F') • 
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Thus 


da' 


3 


16 




3/2 


[- 2(1 - 3 cos 2 I')B2 J + (sin 2 I')B2 2] 


16 


' (1 - e'^Y 


6 (sin 2 I')B 22 ,i + (1 - I ')2 822,2 + (1 + r )2 822 , 3 ] 


5 n'Rg3j^sinI' 


192 a - 2 


(1 - e ' 2 ) 5/2 


[3(1 - 5 cos 2 I')B3 j + (sin 2 I')B3 2] 


7 "' VJ4 


1024 a'3 


'3 (1 ^ e' 2 )’ 


2 (3 - 30 cos^ I ' + 35 cos'* I ') B,^ ^ 


+ 2 sin^I' (i- 7 cos^I')B,I 2 “ 


9 u' J 5 sin I 


61440 a- 4 




[ 30(1 


,,2 - 14cos2 r + 21 cos'* I'jBj , 


+ 5 sin^ T (1 - 9 cos^ I ' ) Bj 2 - 3 (sin'* I ' )Bj 3 ] 

7 ^ ( r 

+ ygg — 1 12 1^- 2(1 - 3 cos^ I') (2 + 3e'^) + 6 sin^I' (2 + 3e'^)cos2h 

+ 30e'^sin^I' cos2g + lSe'^(l-cosI')^ cos (2g- 2h) + 15e’^ (1 ^ cos I')^ cos (2g + 2h)](E'-Z) 


4(1 - 3 cos2 I')Bg j + 6(sin2 I')Bg 2 + 3(1 - cos I')2 Bg 3 + 3(1 + cos I')2 8(2 


( 47 ) 


as 

ai 


2 1 J 2 , , / ^ 

4 77:'o'2i3/2 [-^^ 2.1 ^^ 2 . 2 ) 


(1 


1 "'^/J22 

" 4 ( i - e ' 2 ) 3/2 


in I' [6(C0S I ' )B22 1 + (1 - cos I ' )B22 J “ (1 + cos I' )B22 3_ 


n' R 3 J cos I' 


J_ ■■ J 3 
a' (l-e'^) 


5/2 


[(ll^lScos2l')B3 3 + (sin2l')B3 2] 
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I n' sin I ' cos I ' 


128 3-2 




7/2 


[-10(3 - 7 cos2 I')B^ j + 2(4-7 cos2 I')B^ 2 - (sin^ I')B^ 3] 


1 n' RgS cos I' p 
2048 3,3 (i_ e^^ [2(29 - 126 cos 2 r +105 cos'*I')B3,^ 

+ sin^I' (7-15 cos 2 I')Bj 2 - (sin'* I')Bj jj 

n f 

— p j sinI''|l2[-2cosI' (2 + 3e'^) + 2 cos I ' (2 + 3e ' ^ ) cos 2h 


1 n ' a ' 2 
64 T 


+ lOe cos I ' cos 2g+5e'^ (1-cosI') cos (2g~2h)-5e'^ (l + cosI')cos(2g + 2h)J (E' - I ) 
t 4( cos I ' )Bg, j + 2( cos I ' )Bg 2 + ( 1 - cos I ' )Bg 3 - ( 1 + cos I ' )Bj ’ (48) 


5S2 1 n'R/J 2 

32 (/:%- 2)5/2 


de ' 


2(1 - 3 cos2 I ' ) [l2e' (f ' - !) ^ 3(4 ^ 3e' 2) si 


sinf' '■6e' sin2f' ^-e'^sinSf 




sin^I' [-18e' sin2g-3(4-7e'^)sin(f' t2g) + 3e'^sin(f' - 2g) + 36e' sin(2f' + 2g) 

+ (28 + lle'^) sin(3f' f 2g) * 18e' sin(4f' *2g) *-3e'^sin(5f' /2g)j|- 


^ 2^1 ® ' [34c ' ( f ' - / ) cos 2 h *3(4 • 3 e'-) sin(f' ^ 2 h) *3(4 + 3 e'^) sin(f' - 2 h) 


+ 6e' sin(2f' +2h) ^6e'sin(2f'■■2h)^e'^sin(3f' + 2h)+e'2sin(3f'-2h)j 
t (l-cosI')^[~18e' sin(2g-2h)-3(4-7e'^)sin(f + 2g-2h) + 3e'^sin(f' - 2g+2h) 

+ 36 e' sin(2f' + 2g~2h) + (28 + lle'^)sin(3f' +2g-2h) + 18e' sin(4f' + 2g-2h) 

+ 3e'^sin(5f' +2g~2h)J + (l + cosI')^ [- 18e' sin(2g + 2h)-3(4-7e'^)sin(f' +2g+2h) 

+ 3e'^sin(f' -2g"2h)+36e' sin(2f' +2g + 2h)+(28 + lle'^) sin(3f' +2g + 2h) 

+ 18e' sin(4f' +2g + 2h)+3e'^sin(5f' + 2g + 2h)J^ 
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- I') [48(l + 4e'2) (f' - Z ) sin g + 24(2 + e ' 2 ) cos 2g 

+ 42e' (2 + e' 2) cos (f ' - g) - 72e' (4 + e' 2) cos (f ' + g) + 24e'2 cos (2f ' - g) - 24(3 + 5e' 2) cos (2f ' + g) 

+ 3e'2 cos (3f ' - g) - 2e' (34 + 9e ' 2 ) cos (3f ' + g) - 24e'2 cos (4f ' + g) - 3e' 2 cos (5f ' + g)] 

+ sin2 I' [- 40e'2 cos 3g- 5e'3 cos (f ' - 3g) - lOe' (6 - 5e'2) cos (f ' + 3g) - 40(1 - 3e'2) cos (2f ' + 3g) 

+ 40e' (4 + e' 2) cos (3f ' + 3g) + 20(s + 6e' 2) cos (4f ' + 3g) + 2e' ( 54 + lie ' 2 ) cos ( 5f ' + 3g) 

+ 40e ' 2 cos (6f ' + 3g) + 5e ' 2 cos ( 7 f ' + 3g)J^ 

^2 ■[- 2(3 - 30 cos2 I' + 35 cos'* I' ) [l20e' (4 + 3e' 2 ) ( f ' - / ) + 30(s + 36e ' 2 + 5e '^ ) sin f ' 


n' Re"J4 


2048 a- 2 


(1- e'2)9/2 

+ 240e' (1 t e' 2) sin 2f ' t 5e ' 2 ( 24 + 5e ' 2 ) sin 3f ' + 30e ' 2 sin 4f ' t 3e"' sin 5f '] 

+ 4 sin2 r (1-7 cos2 r ) [l20e' ( 2 + 5e ' ^ ) ( f' - Z ) cos 2g - 20e ' (l4 t 5e ' 2 ) sin 2g + 5e' 2 (48 + 19e ' 2 ) sin ( f ' - 2g) 

+ 10(8 + 144e' 2 + 23e"*) sin ( f ' + 2g) + 200e' (4 t 3e' 2 ) sin ( 2f ' + 2g) + 50e ' 2 sin ( 2f ' - 2g) 

+ 30(8 + 24e'2 + 36"*) sin (3f ' + 2g) * 5e"* sin (3f ' - 2g) - lOe ' ( 34 + 25e ' 2 ) sin (4f ' -t 2g) 

+ e'2(l92^31e'2)sin(5f'+2g) + 50e'2 sin (6f ' + 2g) ^ 5e "* sin (7 f ' * 2g)j 

- sin'* I' [- 350e'2 4g + 35e"* sin (f' -4g)- 35e'2 (24- 13e'2)sin (f' + 4g)- 280e' (4- 5e'2)sin(2f' + 4g) 

- 70(8- 36e'2 - 7e“*)sin (3f'+ 4g) + 700e' (4 + 3e'2) sin (4f' + 4g) + 14(l04 + 252e'2 + 29e“*) sin (5f'+ 4g) 

+ 280e' (8 + 5e' 2) sin (6f ' + 4g) + 5e' 2 (264 + 37e' 2) sin (7f ' + 4g) + 350e'2 sin (8f ' + 4g) 

+ 35e"* sin (9f ' + 


1 


8192 a'2 (1 - e'2)“''2 


{2(1 - 14 cos2 I' + 21 cos'* I') [480(4 + 41e'2 +i8e“*) (f'-Z)sing 
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+ 120 ( 8 + 16e' 2 + 36"*) cos g - 1800e' ( 8 + 12e'2 + e"*) cos (f' + g) + 15e' (336 + 848e'2 + 85e“*) cos (f'- g) 

- 120 (l 6 + 98e'2 + 39e"*) cos (2f' + g) + 660e'2(4 + 3e'^) cos (2f'~ g) - Se' (624 + 1072e'2 + 95e'‘')cos (3f + g) 

+ 40e'3 (23 + 4 e' 2) cos (3f ' - g) - 60e'2 (33 + 21e'2) cos (4f ' + g) + 180e"^ cos (4f ' - g) 

- 24e' 2 (37 + 5 e' 2 j cos (5f' + g) + 15e'® cos (Sf ' - g) - 180e' cos (6f' + g) - 15e'® cos (7f ' + g)j 

+ sin2 I' (1 - 9 cos2 I') 3360e' 2 ( 1 + 2e' 2 ) (f ' - Z ) sin 3g - 140e' 2 (32 + 9e' 2 ) cos 3g 

+ 105e' (I 6 + 208e'2 + 25e“*) cos ( f ' + 3g) - 840e'^ (3 +e'2^ cos (f'- 3g) + 9240e'2 (2 + e'2^cos(2f' +3g) 

- 420e"* cos (2f ' - 3g) + 1400e' (8 + 12e'2 + e"*) cos (3f' + 3g) - 35e'S cos (3f' - 3g) 

+ 840(4 + 19e'2 + 7e"*) cos (4f ' + 3g) + 21e' (304 + 432e'2 + 3Se"') cos (5f' + 3g) 

* 140e' 2 (36 + 17e' 2) cos (6f ' + 3g) + 120e' 2 (l7 + 2 e' 2 ) cos (7f ' + 3g) + 420e ' cos ( 8 f ' + 3g) 

t 35 e '5 cos (9f ' + 3g)] - sin"' I' [- 7566"* cos 5g - 504e'3 ( 5 - 2e ' 2 ) cos ( f ' + 5g) 

- 63e' 2 cos ( f ' - 5g) - 1260e' 2 (4 - 3 e' 2 ) cos (2f ' + Sg) - 105e' (48 - 80e' 2 - I3e' cos (3f ' + 5g) 

504(4 - 25e'2 - ISe"'*) cos (4f ' + 5g) + 1512e' ( 8 + 12e ' 2 t e "> ) cos ( 5f ' + 5g) 

+ 168 (32 + 130e ' 2 + 45e' ) cos (6f ' + 5g) + 45e' ( 240 + 304e' 2 + 23e' cos (7f ' + 5g) 

+ 1260e'2 (7 + 3e'2) cos ( 8 f ' + 5g) + 56e'2 (55 + 7e'2) cos (9f' + 5g) + 756e"' cos (lOf ' + 5g) 

+ 63e ' 2 cos (Ilf ' + 5g)j|- 

1 r r 

^ i I T j |^24e ' |^-2(l-3cos^I')+6sin^I' cos 2h +- 10 sin^ I' cos 2g 

+ 5(1 - cos I' )2 cos (2g - 2h) + 5(1 + cos I')2 cos (2g+ 2 h)J (E' - Z ) 

+ 4(1 - 3 cos2 I' ) 1 ^ 3(4 + 3e'2] sin E' - 6e' sin 2E' + e'2 sin 3E'l 
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I 


+ 2sin^I' j^“9(4 + 3e'^)sin(E'-2h) + 18e ' sin(2E' -2h)-3e'^ sin(3E' "2h) 

- 9(4 + 3e'2) sin (E' + 2h) + 18e' sin (2E' + 2h) - 3e'2 sin (3E' + 2h) 

- 15 1^(2 + 3e '2) - 2 (l-e' 2 ^V 2 +2e'^ sin (E' - 2g) 

+ 6e' |^l-2(l-e'^)’^^ + (l-e'2)“'/^(l + e'2)^sin(2E'-2g) 

- |(2 - 3e'2) - 2(1 - e'2)i/2 + 2e'2 (1 - e' 2 )'i/ 2 j sin (3E' - 2g) 

- 15 -^(2 + 36'^) + 2(l-e'2ji^2 -2e'2 (l-e'^^"^'^2|sin(E' + 2g) 

+ 6e'-^l + 2(l-e'2)i/2 - (i-e'2)"*^^(l + e'2)|sin(2E' +2g) 

- 1(2 - 3e' 2 ) + 2(1 - e' 2 ) 1/2 - 2 e ' 2 (1 - e'2)-i/2| gin (3E' + 2g)j 

’ (l-cosI')2 - 15|(2f3e'2)-2(l-e'2)i/2 + 2e'2 (l-e' 2 )‘i/ 2 ^ 5 jp,^g' -2g+2h) 
+ 6e' [l - 2(1 - e' 2 ) 1/2 + (1 - e' 2 )-i ^'2 + g'2)j sin (2E' - 2g + 2h) 

- '^(2 - 3e'2) - 2(1 - e'2)i/2 + 2e'2 (1 - e'2)‘i/2| sin (3E' - 2g + 2h) 

- 15|(2 + 3e'2) + 2(1 - e'2)i'2 - 2e'2 (l - e ' 2 ) - i'2| gin (E' + 2g ~ 2h) 

+ 6e'|l+2(l-e'2)l/2-|i-e'2)”i/2(i+e'2)J.sin(2E' + 2g“2h) 

- |(2- 3e'2) + 2(1 - e'2)i-'2 - 2e'2 (l - e'2)-l/2| sin (3E' + 2g - 2h)j 

+ (1 + cos I' )2 - 15 |(2 + 3e'2) - 2(1 - e'2)i'2 + 2e'2 (1 - e'2)-l''2| sin (E' - 2g - 2h) 
+ 6e' |l - 2(1 - e'2)i 2 + _ e'2)-i'2 t e'2)j sin (2E' - 2g- 2h) 

- |(2 - 3e'2) - 2(1 - e'2)i''2 f 2e'2 (1 - e'2)"i''2| sin (3E' - 2g - 2h) 

- 15j(2 + 3e'2) + 2(l-e'2)i/2_2e'2^1-e'2)"l'2|,gjp,^£' + 2g+2h) 
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+ 6e' + 2 (l - e' 2 )1/2 _ _ g. 2 ^- 1/2 (1 + e ' 2 )| sin ( 2E' + 2g + 2h) 

- |(2-3e'2) + 2 (1 - e'2)i/2 - 2e'2 (l - e' 2 )‘ i/2| sin (3E' + 2g + 2h) | 


1 n' a2 
384 e L ' 


-yj jl2 [- 2 (1 - 3 cos2 I') (2 + 3e’2) +6 sin2 I' (2 + 3e'2) cos 2h 


+ 30e'2 sin2 I ' cos 2g + 15e'2 (1 - cos I')^ cos (2g- 2h) + lSe'2 (1 + cos I')2 cos (2g + 2h) 

+ 12 ( 1 - 3cos 2 I') jse ' (4 + e ' 2 ) cos E' - 6e ' 2 cos 2E' + e ' ^ cos 3E'J 

+ 18 sin2 I' 1^- 3e' (4 + e'2) cos (E' - 2h) < 6e‘2 cos (2E' - 2h) - e'^ cos (3E' - 2h) 

- 3e' (4 + e' 2) cos (E' + 2h) + 6e' 2 cos (2E' + 2h) - e ' ^ cos ( 3E' + 2h ) 


- 5e' 

{(2 fc'2) 

- 2(1 - e'2)i/2| - 2g) 

*2{( 

2 t e'2) - 

2(l-e'2)i''2 ^lte'2)|cos(2E'-2g) 


[(2-e'2). 

-2(l-e'2)i/2^(,Q3(-3£' -2g) 

- 5e' 

{(2. e'2) 

+ 2(1 - e'2)i''2| cos (E' + 2g) 

+ 2 {( 

2 t e' 2) ^ 

2 (1 - e ' 2 ) 1 / 2 4 c ' 2 )| cos ( 2E' + 2g) 


[(2-c'2) 

t 2 (1 “ e'2 )i/2| cos (3E' ^ 2g)j 

+ 9( 1 

- cos I ' ) 2 

■•5e''^(2 + e'2) “2(l - e'2)i/^J cos(E' “2g + 2h) 

+ 2{( 

;2 + e'2) - 

2(l-e'2)i/2 4 e'2)Jcos(2E' ~2g+ 2h) 


- e' |(2-e'2) - 2(l-e'2)i/2j,(,gg^3£< - 2 g ^ 2 h) 

- 5e' |(2 + e' 2) + 2 (l - e ' 2 ) 1 / 2| ^os (E' + 2g - 2h) 

+ 2|(2+e'2) +2(l-e'2)i/2 + e '2) j cos ( 2E' + 2g - 2h) 
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- e' |(2 - e'2) + 2 (1 - cos (3E' + 2g - 2h)j 

+ 9(1 + cos I ')2 - 5e' |(2 + e'2) - 2(1 - e'2)i^2| (g- _ 2g- 2h) 

+ 2 |(2 +e'2) - 2 (1 - e'2)i/2 (l + e ' cos ( 2E' - 2g - 2h) 

- e' |(2 -e' 2) - 2 (1 - e' cos (3E' - 2g - 2h) 

- Se' |(2 + e'2) +2(1 - e' cos (E' + 2g + 2h) 

+ 2 -^(2+6'^) + 2^1-e'^)^^2 |l + e' 2j^ cos ( 2E' + 2g + 2h) 


“e'|( 2 -e'^) + 2(l"e'^)^^^^cos(3E' ■*- 2 g + 2 h) > — - sin E\ ( 49 ) 


\ 2'^ 3/2 I ' [3e' cos ( f ' ■* 2g) " 3 cos (2f ' + 2g) + e' cos ( 3f ' + 2g)] 

* P 2 j 

+ ^ ~ 1 ' ) ^ COS (f' +2g-2h) + 3 cos ( 2f ’ + 2g - 2h) + e ' cos ( 3f ' + 2g - 2h)] 

+ (l + cosI*)^ [3e 'cos(f' +2g+21i) + 3 cos (2f' + 2g+ 2h) + e' cos (3f' + 2g + 2h)j^ 


J_ n' J 3 sin I 

32 a' (1'- 


(1-5 COS ^ I 


') 12e'(f'-Ocosg + 3e'^sin(f'-g)-*6(2 + e'^) sin(f'+g) 


+ 6e' sin (2f' + g) + e'^ sin (3f' + g)J - sin^ I' j^lSe'^ sin (f' + 3g) + 30e' sin (2f’ + 3g) 
+ 10 (2 + e'2) sin (3f ' + 3g) + ISe' sin (4f ’ + 3g) + 3e ' ^ sin ( 5f ' + 3g)j| 


1 n' ^£"^14 I 


128 „-2 


( 1 -e'^y 


sin I r 

> 2\7 / 2 ■|(l ~ 7 cos^ I ') 60e (f' - 0 sin 2g- lOe'^ cos (f' - 2g) + 30e' (4+ e'^)cos (f' + 2g) 


+ 20 (2 + 3e ' ^) cos ( 2f ' + 2g) + lOe ' (4 + e ' ^} cos ( 3f ' + 2g) + 15e ' ^ cos ( 4f ' + 2g) 2e ' ^ cos ( 5f ' + 2g)J 
~sin2l' j^35e'^cos(f' +4g) + 105e'^cos(2f' + 4g) + 35e '(4 + e'^) cos(3f'+4g) 
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35 (2 + 3e' 2) cos (4f ' + 4g) + 21e' (4 + e'^) cos (Sf ' + 4g) + 35e cos (6f ' + 4g) + 5e' ^ cos (7f ' + 4g)jj- 
-3 _ g, 2 J 9/2 |2(l - 14 cos2 I' + 21 cos'* I') [l20e' (4 + 3e'2)(f ' - i) cos g 


2048 3-3 


+ 30(8 + 24e' ^ + 3e"*) sin ( f ' + g) + 60e' ^ (6 + e' sin ( f ' - g) + 60e' (4 + 3e'2) sin (2f ' + g) 

+ 60e'^ sin (2f ' - g) + 20e'^ (6 + e'^) sin (3f ' + g) + Se"* sin ( 3f ' - g) + 30e ' ^ sin (4f ' + g) 

+ 3e"* sin ( 5f ' + g)J - sin^ I ' (l - 9 cos^ I') [s40e'^ (f ' - 1) cos 3g + 420e' ^ (6 + e ' ^) sin ( f ' + 3g) 
+ lOSe"* sin ( f ' - 3g) + 420e ' (4 + 3e' sin (2f ' + 3g) + 70(g + 24e ' ^+ 30 "*) sin (3f ' + 3g) 

+ 210e ' (4 + 3e' 2 ) sin (4f ' + 3g) + 84e' ^ (6 ^ e' sin (Sf ' + 3g) + 140e' ^ sin (6f ' + 3g) 


■' 15e' “ sin ("7f ' ♦ 3g)J ♦ sin“ I ' |^315e ' si n ( f ' < 5g) + 1260e'^ sin ( 2f ' 1 5g) 

" 420e' 2 (6 + e' sin (3f ' + 5g) + 630e' (4 t 3e' sin (4f ' + 5g) t 126 (8 + 24e' ^ * 3e"*) sin ( 5f ' + 5g) 


* 420e' (4 + 3e'^) sin(6f ' t 5g) + 180e'^(6 * e'^) sin ( 7 f ' + 5g) + 315e ' ^ sin (8f ' ^ 5g) 


35e ' '' sin (9f ' + 5g) 


']} 


1 n ' / ^(2 

64 £ 



-60e'^ j^2sin^I' sin2g + (l - cos I'’)^sin(2g-2h) 


* ( 1 + cos I ' ) ^ sin ( 2g + 2h)J (£'-/)’ 2 sin^ I ' |l5e ' |(2 * e' 

- 2 (1 - e' 2) cos (E' - 2g) - 3 |(2 + e' ^) - 2 (1 - e' (l 4 e ' = )| cos ( 2E' - 2g) 

+ e' |(2-e'2) - 2 (1 - e'2)i/2| ^os (3E' - 2g) 

~15e' -^(2 + 0 '^) + 2(l-0'^)*^^|cos(E'+2g) 

+ 3 |(2 + e' 2 ) + 2 (1 - e' (l + e ' ^ cos ( 2E' + 2g) 

- e' |(2 - e'2) + 2 (1 - 0 ' 2)14 2| cos (3E' + 2g)j 
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+ (l-cosI')^ j\se' |(2 + e'^) - 2^1-e'^)*''^^cos(E'-2g + 2h) 

- 3 |(2 + e'2) - 2 (l - e'2)>/2 (l + e ' =) j cos ( 2E' - 2g + 2h) 

+ e' |(2 -e'2) - 2(1 - e '2 )i/ 2 j (jj-. _ 2 g + 2h) 

- ISe' |(2 +e'=) + 2 (1 cos (E' + 2g - 2h) 

+ 3 |(2 + e'^) + 2 (1 - e' (l + e' cos (2E' + 2g - 2h) 

-e' -^^ 2 - 0 '^] + 2(l-e'^)^'^^|cos(3E'+2g-2h)j 


-* (1 

)■ cos I ' )^ 

|l5e' { 

(2 + e ' 


2 (1 - e ' 

cos (E' - 2g 

-M 


2(1 - e 

'2^1/2 

(1^ 

e ' cos (■ 2E' - 2g - 2h) 

+ e' 

{(2-,'=)- 

■2(1- 

p < 2 ) 1 / 2 ^ cos ( 3E' - 

- 2g - 2h) 

- 15e' |(2 + e'^' 

1 4 2(1 

-e-) 

1 /2 j 

cos (E' 

4 2g 4 2h) 


(2 + e' 4 ; 

2(1 - e 

> 2 y/2 

(1^ 

e'2)| cos (2E' + 2g+ 2h) 


-e'|(2-e'2) + 2(l-e'2)'''2j cos(3E' + 2g+ 2h) 


’ ( 50 ) 


and 


'^^2 _ ^ J 2^ 

^ " 4 2 ) 3/2 


■^Gsin^l' [-2(f'-Z) sin2h-e' cos(f'~2h) + e' cos(f' +2h)] 

- (1- cosl')^ [3e 'cos(f' + 2g-2h) + 3 cos (2f' + 2g“ 2h) + e' cos (3f' + 2g~ 2h)] 

+ ( 1 + cos I ' ) ^ [3e ' cos ( f ' + 2g + 2h) + 3 cos (2f ‘ + 2g + 2h) + e' cos ( 3f ' ^ 2g t- 2h)j^ 

^ -f36 2 sin^ I' (2 + 3e'^) sin 2h + 5e'^(l “ cos I')^ sin (2g“ 2h) 


1 n ' a '^ / 
192 £ 
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”5e'^(l +COS I')^sin(2g + 2h)J (E' ~ i) 


+ 6 sin^ I' |9e' (4 + e'2) cos (E' - 2h) - 9e'^ cos (2E' - 2h) + e ' ^ cos (3E' - 2h) 


- 9e’ (4 + e'2) cos (E' + 2h) + 9e' ^ cos (2E' + 2h) -e'^ cos (3E' + 2h)J 

+ 3( 1 - cos I')2 j^- 15e' |(2 + e'^) - 2 (l - cos (E' - 2g + 2h) 

+ 3 |(2 + e'2) - 2 (l - e' 2)1/2 (l + e ' 2)| cos ( 2E' - 2g + 2h) 

- e' |(2 -e'2) -2 (l -e'2)i/2| ^os (3E' - 2g + 2h) 

+ 15e' |(2 + e' 2) + 2 (l - e ' 2 ) i / = j (•£< + 2g - 2h) 

- 3 |(2 + e'2) + 2 (l - e' 2)1/2 (l + e ' 2 ) j cos ( 2E' t 2g - 2h ) 

+ e' |(2 - e'2) + 2 ^1 - e ' 2 ) 1 ^ 2^ C 3£< + 2g - 2h)J 


3( 1 + cos I' )' 


ISe' |(2 + e'2) - 2 (1 - e'2)i/2j cos (E' - 2g - 2h) 

3 |(2 + e'2) -2(1 -e'2)i/2 + e'2)j ^03 (2E' - 2g - 2h) 

e' |(2 - e'2) - 2 (l - e' 2 )i '2| cos (3E' - 2g - 2h) 

' |(2 + e' 2) ^ 2 (1 - e'2)i"2| (j.. + 2g + 2h) 

|(2 + e' 2) + 2 (l - e' 2)1/2 (l + e ' 2)j cos ( 2E' + 2g + 2h) 

- e' |(2-e'2) + 2(l-e'2)i/2^cos(3E' +2g + 2h)^ 


- 15e 


K 


and 


AL = L - L' . 
AG = G - G' , 

AH = H - H' , 


• ( 51 ) 
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the short-period perturbations in the Keplerian elements a, e, I are obtained by 


a 


e 



+ 



and 


I = r . (52) 

and, in the expressions for AL, AG, AH, the variables / , g, h may be replaced by the variables I ', 
g ' , h' with an error of the fourth order. 

THE LONG-PERIOD TERMS; ELIMINATION OF IT 

At this stage, the Hamiltonian F' - F„' + F/ t depends on the variables g', h' , L' , G', and 
H'; and V is a constant with respect to time. The next step in von Zeipel's method consists of the 
elimination of h' and g'. The elimination of h' is performed by means of a generating function 

S' ■- L' / ' * G' g' H h' - Sj' <• $2 + ••• (53) 

of a canonical transformation from (V , g' , h' , L' , G' , H' ) to (Z" , g' , h' , L " , G" , H"). The new 
Hamiltonian 


F 


•> F, 


(54) 


should be independent of h . 

In order to simplify the formulas, the following symbols are introduced: 


/j - Mg , b Kg , 

T] - — - y 1 - e 2 , 
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and 


e‘ 


ir 


cos I ' . 


Then, from Equation 39, 


F' 


,,2 n„^ L ^ 

/ , * u ' ^ ® 

+ n_ H 


2L'^ 


16/j.^ e 


-^(5- 3 t)' 2 ) j^(- 1 + 3 i 9 ' 2 ) + 3(1 -6I'2) cos 2h'J 


15 (l-f7'^)[-2 + h') + (l-fc''^) cos 2g' + ~2 (l-n>')^cos2(g' - h' )J^ 

T l'6 ^-1 J2 (1 - 3^''^) + 6 b2 Jjj (1 - f^'2) cos 2h'] 


3 1'3 J 

8 L'8'ri“ (1 - - 5H'") sin g' 


3 


128 


i. u j . r n 

[(5- 3t;' 2) (3 - 30fV'2 + 35^^"') - 10(l - (l - 0 '^) (l - 7^^ ' ^ ) cos 2g'J 


lie lM2^ (i-,-^)‘''^[ 6(7-3V^) {l-14«'^^2F<’^)sing' 


7(l-’i'^)(l '• ’ ^ ) ( 1 ■■ 9' ' ' ^ ) sin 3g' j • (55) 

The elimination of h' depends on the solution of the system 


f; (L") 

^ Fo' 

(L-) 

2L"2 


f; (H-) 

= F/ 

(H-) 

= n*H- . 


f; 

^ F.' 

5Sj' 

dh' 

,9H" 

(56) 


and, since S^' does not depend on Z ' , 


L" = L' . 
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As before, a particular solution is given by 


^2 " ^2s “ part of Fj' independent of h ' 


and 


aSj' 

^ ^ {F2 '-f;J . (57) 


Then it is found that 


F 


2 


16/j^ e 


(S- 3 t 7 " 2 ) ( - 1 + 36 »" 2 ) 


+ is(l -27"2 ) (i -0"2) cos 2 g“] - \ b 2 ^ (l - 3 f ^"2 ) 

“ I J3 r .8 „.-5 ■ 2,- 2)1^2 (1 -r"2)>''2 (l - 5,.' 2j gin g 

L, T] 

- m b' J4 L- 10'^. 7 [( 5 - 3 ^' 2 )( 3 - 30 ('-^+ 35 H'")- 10 (l-r,'' 2 ) (j - ^ 2 ) ( i _ y,." 2 ) ^os 2 g"] 


5 

256 




79 U 


( 1 - 77 ' 2 ) 1'2 2 ) 1/2 [g (7 .. 3 t/ 2 )(i- 14 - -2 4 21 f- 


sin g 


7 (1 - 77 '^ ) ( 1 ^ ^ ) ( 1 ^ 9- 2 ) sin 3g' 


(58) 


where 


v" 


L' 


and 


0 


H' 


G" 
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Furthermore 


1 f 3 / ,, , 

( 5 - 377 " 2 ) (l-e" 2 ) Sin2h 


+ ^ (l-^"0 [(1 sin2(g' +h')-(l-6")2 sin2(g' -h')]} 

- -|n'2b2 ^ (l -0"2)sin2h' . 

The long-period perturbations depending on the motion of the node are given by 


L" + ^ L“ 

ai ' 




\ 1 - 1 





■ (1 ' 

r " cos 2(g ' * h" ) - (1-6 

i COS 2 (g" - h )] ^ 

/ 



tIt [(s- 3 r, ' 2 )(i ■ ,^" 2 ) cos 2 h- (l 77"^ ) {(1 ♦ cos 2 (g" + h") 

+ (1 ---)^ cos2(g" -h' )}] l^^cos2h''j 


H 1 2h" + 15(1 - ^ v“^) [(1 + 0")' si 


sin 2( g ' + h" ) 




'll 9 n ' ^ 

!in2(g" -h")]^ - -4 — Jj 
J 


1 - 5"2 


sin 2h" , 
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5G" 


^ e" 16 e 


^ [3(50"2 - 3t?'' 2) sin 2h" - ^ (1 +d"){v“^ + ^") sin 2(g" +h") 


+ ^ (1 -0") (t 7"2 -0") sin2(g" -h")] - J 22 


3 - 50" 


sin 2h" 


( 64 ) 


and 


5H" 


h" + (3'5" (3-^"^ - 5) sin2h" + ^ ( 1 - 77 " ^ ) [ ( 1 + 0 " ) sin 2(g" + h" ) 

+ (1-0") sin 2(g" - h")]j 


3 n ^ r . , 

2 — ^ J22 T - -'4 2h 

n. L V 


( 65 ) 


It is important to note that these terms are factored by first-order factors. Then, a second- 
order theory produces first-order long-period terms. It is necessary to go to third order to ob- 
tain second-order long-period terms. This is done next. 

THE SECOND-ORDER LONG-PERIOD TERMS 
AND ELIMINATION OF AND THE TIME 

The small parameters of third order are as follows: 


‘(£ 
2 I n 


(earth's oblateness) 


(sun's perturbation) 


n^y i(j (correction due to the inclination of 
sin , 

the moon's orbit to its equator) 


(solar radiation) 
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and 


(eccentricity of the moon’s orbit) 

(earth's perturbations— third 
Legendre polynomial) 


an 

Q 


n 


(correction due to physical libration) 


Since short-period perturbations of the third order will be neglected, it is understood that all 
terms of third order which depend on l will be ignored. Suppose that is the third-order part of 
the Hamiltonian and F^' the part free of short-period terms. The Hamiltonian will be time- 
dependent through the longitudes of the earth and of the sun. This fact introduces one more degree 
of freedom in the problem, and it is necessary to introduce the time as a canonical variable. This 
is done through use of the following equations { r - t ^ const.): 

d; d'‘> 

Ht ^ ■ FT ■ 1 (66) 


and 


dt 


£3 ^ ^ 

dr - dt ' dt ' 


( 67 ) 


so that 


3 F - T 


(68) 


and 


T - Fj + const. 


( 69 ) 


This additional complication does not affect the development of the theory up to this point, if 
one sets 



( 70 ) 


The next step is then to eliminate h' to third order to obtain second-order long-period per- 
turbations depending on the argument h' . 
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The orders of magnitude are as follows: 


3o' = F ; - T = — - T 


3i' = Fi' = "eH' 


3a = 


and 


(71) 


where T has been incorporated into 3^' . The von Zeipel differential equation of the third order is 
then 


or 


3s; aSj' aSj' as^ a'3^ 

3 ah' aH" dg' ac" ah' aH" a-r ax 


aSj' as/ 

a G“ a g ' 


as; 

aF/ 

aF/ 

as; 

* ^ 


(9h' 

aH" 

aG" 

ag' 

ah' 

dr 


aF/ as/ 

F + — - 

^ dg' ac" 


(72) 


(73) 


Since we had supposed thatas/Zar = 0 in the previous order evaluation, a particular solution 
of the third-order equation is found by taking 


F 


3 



where the subscript s designates non-dependence upon h" andr. Further, 


and 


as, 

ah 


- F 


3ph 


aH"4^. 


,aG" ag' /pj,' \ag' aG" /p^/ 



F ' 

3pr 


as/ 





\aG" 





(74) 


(75) 


( 76 ) 
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where the subscript ph' indicates the inclusion of all trigonometric terms with arguments of the form 
jh' + kg', where j and k are integers (j / 0), and the subscript pr indicates the inclusion of all trigo- 
nometric terms with arguments of the form jr + kh' + mg' , where j , k, and m are integers ( j / 0). 

Then S^' will be obtained as a periodic fimction independent of the time, and S3' as a periodic 
function dependent on time. 

Since the functions F^'p, F^', S^', and F^" are known, the only parts that have to be computed 
are those corresponding to the third-order terms of the Hamiltonian (F3' ) . 


THIRD-ORDER TERMS GENERATED BY COUPLING OF SECOND-ORDER TERMS 

It was found (from Equations 55 and 58) that 


^ 2 * 16 


§7 n'L' [^(- 1 + 36 ^'^) ^3(1-6''^) cos 2 h'] 


15 


(l-’)'^)['2 (1 ' *^')^cos 2 (g' th') + (l-«' 2 ) cos 2 g' + -j (1 - 6 '')^ cos 2 ( g' ' )]|- 


■? L -6 -3 [■ J2 (1 ■ 3 f'' 2 )* 6 b 2 J32 (1 - f-'2)cos 2 h']-| — i— T (1 (1 _ 1 


b'* J 


128 L' 10 


77 [(5 - 3 ^' 2) (3 - 30 i -‘^ t 35 f^"*) - 10(1 - r,' 2) (1 - «'2) (1 - 7f^'2) cos 2g'] 


5 

256 L' 12 9 


(1- (i-.',-'2^i/2 ^ 2 P/'’) sin g' - 7(1- 7'^) (l-;.<'2)(l-9t?'2) sin 3 g' , (77) 


16? (■^) L' [(5 - 377' (- 1 + 3 r '^) + 15(1 - v "^) (1 - cos 2g'J 


^ b2 ^ (1 - 3ft='-2) - I b3 J3 g (1 -77 ''2) 1/2 (1 ^ 6>"2 )i/ 2 (1- st 7 ~ 2 j sin g" 

3 ^ 

1 _ 4 -r 

128 


J4 ,^,10 ,.7" [(5 " 37" 2 ) (3- 3 O 0-2 + 356 '"'*) - 10(1-7"^) (1 - (1 - 76 i" 2 ) nos 2 g"] 


sis (i-7"^)'"2 (1-0"2)V2 [5(7. 3^ ..2) (1- 140"2 + 210"") sing" 


- 7 1 - 


(1 - 7"^) ( 1 - 0 "^) (1 - 90 " 2 ) sin 3 g"j ’ ( 78 ) 
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and n* Sj' is given by Equation 59. In the coupling terms, the variables in F^' can be double-primed 
with an error that is of the fourth order. 

The partials needed are given below. 


dH" 


3 

16e " 


- 7 ^ ( 2 ( 5 - 3 t]" 2 ) 0"2 (1 - cos 2h") 

16 ^ 


+ 5(1 - T)"^)e'' [(1 + 6") cos 2(g" + h") - 2B" cos 2g" 

1 3 

- (l-e' )cos2(g"-h")]} + 2 b2n'2-^-^ (J^ ~ 2 J 2 ^ cos2h") 


3 , ^n' ^ f?" (11 - 15f^"2) (1 - 

8 ^ J 3 (1 - ^" 2 ) 1/2 


sin g" 


15 


pt2 n ' 2 


32 J 4 L' 5 


[( 5 - 3 ^ 2 ^ ( 3 _ 7 ,,- 2 ) - 2(1 - r,' 2 ) (4 - 7f'"2) cos 2 g ] 


15 

256 


-■Js [ 2 ( 7 -3,-.) (2,- .2.-7..03-. •) 


sin g 


- 7(1 - • 4 " 2 ) (1 - f^"2)(y - 15fr"2) sin3g"j ’ ( 79 ) 


dC 


3n 

TFi 


- (- 7 ) "T [ 2 ( 7 "^ - 5 r 2 ) + 2 ( 5 " 2 - 3 t)" 2 ) cos 2 h' 


5(1 + t-' ) ( 7"2 ^0‘) cos 2(g" + h' ) + 10 (:-' 2 - t"'' 2 ) cos 2g" 

5 ( 1 - 0 ") ( 0 " - V"^) cos 2 (g" -h")] + 1-77^ [jj (1- 50"2) - (3- 50 " 2 ) cos 2 h"] 


I b2j 




3 L'3t7"6 (1 -^”2)1/2 (i_e..2)V2 


[ 5 - 410”2 + 400'"’ - 7 ,"^ (4 - 350 "2 + 35e"") sin g" 


b" I 2 I - 2700" 2 + SSSf-'"* - 9v"^ (1 - 14f-' 2 + 21- '*) 


- 2 [y - 720"2 + 770'"' - 77 " 2 (s - 560" 2 + 630'"')] cos 2g"| 
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+ S b'Js {2[7(9-164a"2+4410"^- 2940"«) 

- t?"2 (77 - I445t?"2 + agSBe'"* - 26676'"®) + 377'"* (6- 1196I"2 + 3360"“ - 2310"®)] sin g" 

- 7(l-r,"2) (l-0"2) [3-370-2 +420 "“-t 7"2 (2- 270-2 + 330"“)] sin 3g"} 




sin 2g" 


I 3^3^ /Tn'2 (i-^" 2)V2 (i_^..2)V2 (i_5y..2) 

8 J 3 r.2„..5 


- U4 T n'2 (1 - ,,"2) (1 ^ .r 2) (1 _ y,.-2) 

32 J4 T '4 „..7 sin 2g 


•I . c 

256 '2® J, 


l'2 (1 - n'2)l'2 (1 - , 2)1 


[2(7 - 3t"2) (29 126" 2 . 105 - '“) COS g 


7 (i ->7 '2 )(i . 2) (7 - 15,. 2) 3 g •] 


0h' n* ^ 2 P 


^ -2)<,os2h 


7 5 (1 n 2 ) [(1 + .. )2 cos 2(g 7 h ) * (1 - f )2 cos 2(g - h )]| 


3 n ' 2 1 - 77 " 2 

2 ° J 7— cos 2h ' 

"c 0 ^ 


15 / \ r 

j (1-77"2)L' [(1+77")2 cos 2(g" + h ' ) - ( 1 - ti" )2 cos 2(g" - h" )] 


32® Vn^ 


3^ (^^') ^ - 377-2) sin 2h" - 5(l + 0")(r,"2 + ^.. ) sin2(g" + h" ) 
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+ 5(1-6'") {ri" 2 - 6 " ) sin 2( g" ~ h" ) j “ 4 




(- 3 + 56" ^ ) sin 2h" 


(84) 


From these equations and Equations 74 and 75, it can be shown that 


Sj' (coupling) 



326" (1 77-2 


(15 ' 1777" 2 ) sin 2h 


- At -" (1 - 6"2) (5 - 377” 2) 2 sin Ah " 

+ 80(1 + ;+")2 (2 - 36")t7"2 (1 - t 7"2) sin (2g' + 2h") 


+ 80(1 - f >")- (2 36 )77 2 (1 - t;" 2^ sin (2g" - 2h") 


+ 10(1 + 6”)2 (1-6") (1 - 76 2^ [^5(1 - - 57;" 2 j sin (2g ' + 4h" ) 

+ 10(1 - (1 7 ) (1 - +1 2) [S(l t -6 ) - 67 , 2] sin (2g - 4h ) 

+ 25(1 7 6" )2 (1 - ti"2) [(1 - - 27.,"2] sin (4g 7 4h") 

+ 25(1 - - ')2 (1 - 77 2) 7 f.' ) - 277 '2j sin (4g' - 4h )| 


9 

256e 



2 ) (5 - 3> ' 2) sin 2h 


+ 5(1 + 6")2 (1 7 2 ^-" - 56"2) (1 - 77’ 2) sin (2g' + 2h' ) 

+ 5(1 - 6")2 (1 - 2+’ - 5-' 2) (1 - 77" 2 ) sin (2g ~ 2h" ) j 


9 

512e 


J,,b2 {-8- (1 


“"2 ) (5 - 37- 2) '2 sin 2h" - sin 4h 


- 20(1 + 6")2 (1 - 6 ") (3 - 5‘ ') (1 - ti- 2] sin (2g" ^ 2h") 


- 20(1 - 6")2 (1 + 6' ) (3 + 56") (l - 77’ 2) sin (2g" - 2h") 
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+ 5(1+0 ")^(l“(9")(l + 56'") sin(2g"+4h") 


5(1 - 0")^ (1 + 0") (1 -50") (l - v'^) sin (2g" - 4h")^ 


9 /n'V .0"(l-0"^) 9 /n'V 0"(l-0"2) 

^ J. J. ' - 16 y 4^ JA sin4h" 


9 /xn' (i- 7]"2)V2 (i_ ^..2)1/2 /n^Y r ^ 

" J3bM(l+«") [25(1-0" -0"2- 70-3) 


- t 7"2 (9 + 150" + 150"2 - 1750"^)] cos (g" + 2 h") 

+ ( 1 - - ') [25 (l + - - '2 + 7 ,. 3) 

- 7 ^ (9 - 15(> ' t 15; - ' ^ + 175' - '•’ )j cos ( g" ■ 2 h ' ) 

- 25(1 - 7 ' 2 ) (1 • •• )2 (1 . 2 (>- • 7 "' 2 ) cos ( 3 g' t 2 h") 

- 25 (1 - n" ^ ) ( 1 - 0 " (1 - 20" - 70 " ^ ) cos ( 3g" - 2h" )!■ 

9 / n ' V un ' ^ (1 - Tj " 2 ^ 1 ' 2 ^ _ fy 2 \ 1/2 r , . 

" T28 t '3 9 J 3 J 22 b= (1 * 0 ') ( 3 - 250" + 50-2 +250'2) cos(g" + 2 h") 

\"£/ V 

^ (1 •" ) (3 - 25 - + 5 -"-^ - 25 c 2 ) cos (g - 2 h")] 

J4 {s'' (1 " ■' [35 - 107'2 (2 t 7 f 2) - 3 r,'+ (1 ^ 140 2 )j sin 2 h" 


45 n' 2 n\ 
81'92. L- ^ 


+ (1 - 7 2) + „-)2 [35(3 + 4f, ■ - 3O0"2 - 200"^ + 510'"* ) 

- 3t 7"2 (11 + 360" - 1260"2 - 1960"^ +3150'"*)] sin (2g" + 2h" ) 

+ (1 - t7"2) (1 - 0")2 [35(3 - 40" - 300 "2 + 2O0"2 + 510'"* ) 

- 3t]" 2 (11 - 360" - 1260"2 + 1960"3 + 3150'"*)] sin (2g" - 2h") 
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- 35(1-77"2)2 j^(l+0")2 (l + 26?" - 90"2) sin(4g“ + 2h") 

+ (l-i9")2 - 20 " - 90 " sin (4g“ - 2h" 

45 /n ' y n ' ^ ( 1 - g " ^ ) 


512 I n* 

^ E/ 


L'2 7,"11 


J4J22b« [4(5-3 t,"2) 0 " { 3 - 70 "^) Sin2h" 


(l - t 7"2 j (1 + 0") (3_ 19^" _ 7^.-2 + 350"3j sin (2g" + 2h") 


+ (l - t 7"2] (1 - r^") (3 + 196?" - 70 "^ - 356?" 3) sin (2g" - 2h")j 


45 n-3 (i - 77 -2) V2 (i _ c^" 2)V2 /n, 


32768e 




J 5 bs|(l+ 6 ?") [35(15 + 416 ^" - 26 26 ?" = 


- 3466?"3 + 5196?'"' t 321r"5 ) + 14 Y ^ (259 - 2556?" - 5666?"=^ + 17106?"^ - 52 56^'"' - 14556?' =) 

- 77'"' (1559 - 10156"" - 54466;"^ + 7350t? + 7357?'"' - 5775‘7"S)J cos (g" + 2h" ) 

+ (1-6?") [35(15 - 416?" - 26 26?"^ -t 3466?"3 + 5196?'"* - 3216?"5) 

6 1477"2 (259 + 2556?" - 5666>"2 - 17100 "^ - 5256?'"' + 14556>"5) 

- 7,'"' (1559 + 10156?" - 5446fl"^-73506? "3 + 735??'"' + 57756?"3)j cos (g" - 2h") 

- 14(i+0"3) [s( 9 - 76?" - 1196?" 3 + 156?" 3 4 174/6'"') - •n"2 ( 179 _ 14^7" _ i3206?"2 + 30i9"3 4 16056?'"') 

+ 3 v '“' (I8+ 70 " - 1356?"2 - 156?"3 4 1656''"')] cos (3g' + 2h") 

- 14(1- 6")3 [5(9 4 70 " - 1190"3 - 156'"3 4 1746''"') 

- 77"^ (179 + 146" - 13206"2 - 306"3 4 16056'"') 

+ 377'"' (18 - 76" - 13 56? "3 4 15/?"3 + 1656"^)] cos (3g" - 2h") 

(1-7,"2) (i-0"2) [i-96"2-t7"2 (3- 116"2)] 

[(1 + 6")3 cos (5g" + 2h") + (1 - 6")2 cos (5g" - 2h")]j. 


+ 105 
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45 

4096 


4Y 




1/2 


^SL‘ v"' 


Js J 22 b’|2(7-3T7"2) ( 29 - 126t^ 2 + lOS "*) [( 


(l + f/'')(3-St^ ")cos (g"+ 2h' ) 


(1 - (9") (3 + 5&") cos (g" - 2h")j - 7 (1 - t)" 2) (l-e/"2) (7-156'"2) j^( 1 + " ) (3 - 50" ) cos ( 3g" + 2h" ) 

+ (l-0")(3+50")cos(3g" -2h")]| ’ (85) 


Q ( 1^(2^ f ^(2 \ f” 

Fg" (coupling) = 1 — r) ( — *"]“ Vl - e"^ cos I " I (2 + 33e" 


32 


- (2 - 17e" cos^’I" + 15e" ^ sin^ I" cos 2g"] 

/ n ' 2 sin^ I " cos I ' T / o\ ,, o ,1 

[2(243e 2) il5e = cos 2g ] 


Q/n'Xn'^sin^I cosi 


4 * . 

\n^/ a 


— T 2 b'* 
- j.- 2^/2 J 22 


( 86 ) 


and, of course, 


S ' (coupl ing) 0 


( 87 ) 


It is important to note that S^' (coupling) is made up of second-order terms, although it is ob- 
tained by multiplying a second-order quantity by a first-order quantity. 

The next few sections are devoted to the computation of F3' . From that computation various 
additional parts of F3" , S3' , and S^' are then derived. 


THE RADIATION PRESSURE 

For this computation, the shadow effect will be neglected. Its inclusion would require, of 
course, the introduction of a shadow function due both to the moon and to the earth. 

Consider o to be the absolute value of the acceleration of the orbiter arising from the solar 
radiation. Then, the disturbing function for the radiation pressure will be given by (see the earlier 
section, "The Angles s^'^ and 8^3 ") 


R - - a V cos S,", , 

O 13 

or, neglecting the inclination of the lunar equator to the sun's orbit, 

= -or fcos ( f + g) cos (n - - sin ( f + g) sin (O - cos I 


(88) 


( 89 ) 
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The elimination of short-period terms gives 


r_ 

2 


cos (g' + h' + A-g, - + sin^ 


r_ 

2 


(e' 


+ A. “ 

0 € 


(90) 


where the node n was written in terms of the canonical variable h. There is no secular contribu- 
tion from this effect, that is, F 3 ; (radiation) = = 0. On the other hand, since there are no 

terms strictly independent of time, the contribution to S^' is zero, or 


$ 2 * (radiation) - 0 - 


(91) 


Therefore, 


5S3' ( radiation) 
dr 


^3pt ®o' 


or 


S3' (radiation) 


2 

2 


o a" e" 



I 

^2 sin 


(e 


■t h' 


- A-,) 

o ) 


I" 

2 sin 


(g' - h' 



(92) 


where 


\ - n T + const. 

0 e 


and n* and n* are, respectively, the mean motion in longitude of the sun around the moon and of 
the moon around the earth. Of course, n* - n* :>r n* is a good approximation. 


THE SECOND LEGENDRE POLYNOMIAL FOR THE 
SUN'S GRAVITATIONAL PERTURBATIONS 

The most important term in the disturbing function due to the sun's gravity is 

P 3 (cos Sj' 3 ) , (see Equation 20) 
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where, for the computation of S^" , the sun can be considered to move along the moon's equator. 
Thus, 

cos Sj'j = cos ( f + g) cos (n - - sin ( f + g) sin (O - cos I . 

On the other hand, neglecting the mass of the moon. 


^ a 3 - k2 (m^ + m^) - (l 


and since the ratio can be neglected. 


Furthermore, if the eccentricity of the earth orbit is neglected, 

FjCsun) - n^2 r2p^ (cos Sj") 

~2 n^2 r2 (3 cos 2 - l) 

- ~2 ^cos2 ( f + g) cos2 (n - - 2 sin ( f + g) cos ( f + g) sin (0 - cos (j - cos I 

t sin2 (f + g) sin2 (q- cos2 l] - y n^2 r2 . 

Rewritten, Equation 93 becomes 


F (sun) 


— n 2 r2 
2 "e 


(cos2 f cos2 g + sin2 f sin2 cos2 (U - 


- 2 (“ sin2 f sin g cos g + cos2 f sin g cos g^ sin (O - cos (O - A^] cos I 

1 1 3 

+ (sin2 f cos2 g + cos2 f sin2 g^ sin2 (Q - A^J cos2 I - ^ ~2 sin f cos f ■ Q , 

where Q is independent of f. It follows that 

Fj'(sun) = n^2 a"2|l 2(2 + 3e" 2 ) (l - 3 cos2 I ") + 30e" 2 sin2 I " cos 2g" 

+ 6 (2 + 3e" 2 ) sin2 I " cos 2 (h" + “ ^ 0 ) + 15e" 2 ( 1 + COS I" cos 2 (g" + h" ^0 “ '^q) 


(94) 


+ ISe " 2 ( 1 - cos I" )2 cos 2 (g" - h" ~ ^ ^o 


.)] 


(95) 
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Ill 


The secular contribution is 

Fj" (sun) = ^2 + 3e"2) (l - 3 cos^ I") + 15e" ^ (l - cos^ I") cos 2g"] • (g0) 

There is no contribution to S^' , and the contribution to S 3 ' is 

S3' (sun) = 64 [2(2 +3e"2) sin2 I" sin 2(h" 

+ Se"2 (1 + cos I')2 sin 2(h" + g" -^g) * 5e'’ 2 ( 1 - cos I" )= sin 2 (h" - g" + \^ - \^)] • (97) 


THE THIRD LEGENDRE POLYNOMIAL FOR THE 
EARTH’S GRAVITATIONAL PERTURBATIONS 

From Equation 20, 

F3 (earth) - {^) P3 ^I'o) ' ( 98 ) 


Now, 


n 


2 

<z 


('^e + "’e) 



^2 


£ 


SO that, considering the moon’s orbit to be circular, 


"c r3 


F3 (earth) = ^7^7^ P3 ' it ^3 ' 2 (S S^;, ' 3 cos S^;,) . 


If 


and 


•^ [( 1 + cos I ) cos (g + h) + (1“ cos I ) cos ( g - h)] 


[(1 + cos I) sin (g + h) + (1 - cos I) sin (g - h)] , 


then 


cos cos f - sin f 
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so that 


F, ( earth) - "o 




3A„^ B- cos^ f sin f 


+ f f " Bgf sin^ f) - 3(A^ cos f - B^, sin f)] • ( 100 ) 


It follows that 


F 3 ' (earth) ~ ~ ^ 


{s [(3 + 4e ' 2) a; 3 + 3A; B^' ^ (l - e ' ^)] - 3A; (4 t 3e ' =)} . (101) 


where 

Aq,'^ ” ^ Ip ( 1 - 0 '^) (1-/9') COS (g' -3h') t(l-0')^ cos (3g' - 3h ' ) 

+ 3 (1 - 0'2) r 1 - 0' ) cos (3g' - h' ) 4 3( 1 - 0') (3 t 20' t 30'^) cos (g' - h') 

+ 3(1 + (!') ( 3 - 2 /'' t 30' =) cos (g' ' h') I 3 (1 -0'2) (1 e 0') cos (3g' + h') 

^ 3 (1-0' 2 ) (1 +0') cos(g' * 3h') + (1 +0')3 cos (3g' ' 3h')J , 

A* ‘ y (1 t 0') cos(g' +h') t I (1 -0') cos (g' -h') , 

and 

A^.'b;2 ^ [(1 * 0') ( 3 - 20' ^ 30'2) cos (g' th') 

t (1 - 0' ) (3 I 20' t 30' 2 ) cos (g' -h') - 3(1 -0') (1 - 0'2) cos (3g' -h') 

+ ( 1 - 0') ( 1 - 0 '^) cos (g' - 3h') - (1 - 0')3 cos (3g' - 3h') 

- ( 1 t 0')3 cos (3g' + 3h') - 3( 1 i O') (1-0'^) cos (3g' + h') 

+ (1 + 0') (1 -0'2) cos (g' + 3h')] . 

The contributions to F/ and to S3' are both zero. The contribution to s^' is then given by 



( earth) dh ' , 


Sj' ( earth) 
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y 


or 


// 3 « 

S^'Cearth) = ft) / [-9(l + cosI") (l + 10 cos I " - 15 cos ^ I ") (4 + 3e"^ ) sin ( g“ + h") 

\"e / ® ^ 

+ 9(l-cosI") (l - 10 cos I" - IScos^ I") (4 + 3e"^) sin (g" -h") 

+ 15 sin^ I" (1 + cos I") (4 + 3e"^) sin (g" +3h") - 15 sin^ I" ( 1 - cos I") (4 + 3e"^) sin (g" - 3h" 
+ 315 sin^ I" ( 1 + cos I") e"^ sin (3g" + h") - 315 sin^ I " ( 1 - cos I " ) e" ^ sin ( 3g" - h" ) 

+ 35( 1 + cos I")® e"^ sin (3g" + 3h") - 35( 1 - cos I ")^ e"^ sin (3g" - 3h")J • (102) 


THE ECCENTRICITY OF THE MOON’S ORBIT 

The correction for the eccentricity of the moon's orbit is given by 


^3 (^c) 


or, keeping only the first power of e^, 


- 1 


(cosSjo) 


^3 (®c) " T e^-f- r^cos (cosSio) 


(103) 


3 “e o ~ 

J e (3cos^ Sj^ - 1) cos Zg 


where the mean anomaly I ^ is I g, ^ t + const. Then, 


Fa' (^c) 


3 ^ 

64 e 


(104) 


2(2 + 3e' ^) (l - 3 cos^ 1 ') + 30e ' ^ sin^ I' cos 2g' + 6 (2 + 3e' sin^ I ' cos 


2 h' 


+ 15e' ^ ( 1 + cos I' cos (2g' + 2h') + 15e'^ (1 - cos I ' )^ cos (2g' - 2h')J cos Zg, . (105) 

There are no contributions to S^' or F/ . The contribution to S3' is 

// 2 

^3 (®e) ^ ^ - ..g — - j^-2(2+3e"^) (1 - 3 cos^ I") + 30e"^ sin^ I" cos 2g" + 6 (2 + 3e"^) sin^ I" cos 2h" 

+ 15e"^ (1 + cos I")^ cos (2g" + 2h") + 15e"^ (1 - cos I")^ cos (2g" - 2h")j sin Z^.(106) 
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It is worthwhile to note that, although this result is a part of a third-order generating function, 
its actual order is e^, , which is considered to be a second-order quantity, because a small 

divisor is introduced through the integration. 

THE INCLINATION OF THE MOON’S ORBIT TO ITS EQUATOR 

The expression for the earth’s perturbation is 

— Pj (cos S,'o) , 


where cos s/g = K(ig). By a Taylor expansion, if sin is small. 


cos S/o - K(0) + 


1 1 

<y 

i 


^ cos S/g 

S (sin i^^ 

sin i^ - cos Sj'g ^ 

cos ig(5i(j 


(107) 




Thus the correction to be introduced is 


^3 (’c) 


<') cos Sj'q 


sin i g CO s S j g 


J ‘c" 


It is easily shown that the bracketed quantity in Equation 108 is 


(108) 


sini sinv.sin(f ■> g) . 


Thus, 


^3 (’e) ' ^ g— cos S,'g sin ( f t g) sin sin I sin ij 


(109) 


Now, 


cosSj'Qsin(f + g) - sin(f + g)cos(f + g)cosh “ sin^(f + g)sinhcosI . 


Therefore, the part of F 3 (i^j,) which is free from short-period terms is 


^3 (’ffi) 


8e 


jin Vjj sin I' sin i^j 2 (2+ 3e'^) cos I' sinh' 


Se ' ^ ( 1 + cos I ' ) sin ( 2g' + h ' ) + 5e ' ^ ( 1 - COS I ' ) sin ( 2g ' - h ' )j . (110) 
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There is no contribution to F^" and S^' . The contribution to S3' is 


S3 


) 


8 c (rig, +N,^ 


sin I' sin - 2 (2 + 3e"^) cos I" sin h" 


+ 5e"^ (1 + cos I") sin (2g" + h") + Se"^ (1 - cos I") sin (2g" - h")j , (111) 


where N„ is given by t + const . Again, this is a second-order contribution, since 


THE NON-SPHERICITY OF THE POTENTIAL FIELD OF THE EARTH 

Because of the fact that the earth's equator is not the reference plane, the form of the dis- 
turbing function due to the zonal harmonic coefficient j ^ of the earth will be derived from basic 
relations. 

In Figure 5, E is the earth, M the moon, and S the orbiter. The disturbing function for the 
motion of s is given by 



m. 


F. 


2 

e J 2 

1 7 ^ ?2 ( sin cp) 




( 112 ) 




Figure 5-Earth -moon-orblt-er configuration. 


where is the equatorial radius of the earth 
and (p the latitude of the orbiter with respect 
to the earth's equator. The part has 

already been taken into account, and the rest 
is a third-order quantity. Therefore, 


k2, 


F3(»-) 


TT R/ 3 ■ 


(113) 


If the terms and are neglected, it follows that 


F,(^) 


"c Re^i2 


P (sin <p) 


(114) 


assuming ^ = 1. The angle ^ must now be expressed in terms of the orbital elements referred 
to the lunar equator. 

In Figure 6 the geometry of the problem is given. It follows that 


sin (p = ( sin I cos 1^ - cos ip cos I sin sin ( f + g) - sin 1^^ sin cos ( f + g) 
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Then, if 


A - sin I cos Ijg - cos 0 cos I sin 

and 


B - - sin sin 0 , 


we have 



Figure 6 — Planar configurations. 


F (®) 


"c i2 T3 


[ 2 13 3 "] 

■4 ( A^ + B^) - 2^4 (b^ - A^) cos 2 (f+g) + -2ABsin2(f + g)J . 


Using the relation 


d/ = df -( 


2 ^ (f * /l - 


cos kf( -1)** df 


where 


i (t -yi -e^) 


(115) 


elimination of short-period terms from Equation 115 yields 


F ' (1.) 


1 2 u 2 : 




A2 t B^ 


/J' ^ ^ - e' 2 ^ ( B^ - A^) cos 2 g' + 2 AB sin 2 g‘ 


(116) 


where 


3 ' " /i(e') . 

Furthermore, since 0 h + \^ -n' , where w' is the longitude of the descending node of the moon's 
equator on the earth's equator, we have 

B^ + A^ = ^sin^ I' cos^ ^ cos^ I' sin^ ^ ^ 

- -2 sin^ sin^ I ' cos 2(h' + Aj, ~ H') 

“ -2 sin2I^sin2I' cos(h' +A^-n') 
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= (-^ sin^ Ig, - sin^ I' cos^ - -5 cos^ I' sin^ I^] 

- -j sin^ Iffi (1 ^' ) ■'■ ) 

+ Y sin 21 ' sin 2 I|j cos (h' + “ ii‘ ) 


and 


■*■ O / • 

AB = 2 ' sin'' I^cos I sin 


2(h' + \^ -n') - ^ sin 21^ sin I ' sin (h' + - O') 


Thus the contribution to S ' is zero, and 


1 "e i 2 K,:* " 

8 t 


F,' (e) 


- (l-3cos^I') 

+ 3/i"^ (l + 2^1 - e"^) sin^ I" cos 2g" ( 1 - 3 cos^ I . (117) 


The contribution to s,' will be given by the integration of 


5S ' (e) 3 n ^2 r ^2 


dr 


16 


2 sin I" [sin I" sin cos 2(h' + ‘ SI') * 4 cos I" cos cos(h' + - O' 


' p' ^ [1 + 2^1 - e" ^ )[- ( 1 + cos I sin I,,, cos 2(h' + g' t - fi') 

- (1 - cos I")^ sin cos 2(h' - g' ^ - ii') 1 4 sin I" (1 + cos I") cos cos (h' + 2g' + Xj. - 0' ) 


- 4 sin I' ( 1 - cos I " ) cos cos 


(h' - 2 g' +x^-n')]| • ( 118 ) 


Since 






we have 


S 3 ' (®) ' 32 




2 sin I" [sin I" sin 1^, sin 2(h" + Xg, - H') + 8 cos I '' cos sin (h" + X^ - 0')j 
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+ /3" ^ ^1 + 2 i\- e" J^-(l + cosl")^ sinl^sin2(h" + 

~ (1“ cosl")^ sinig sin2(h" - g" +8 sini" (1 + cosi") cosl^ sin (h" +2g" 

- 8 sin I" (1 - cos I") cos sin (h" - 2g" + “ ^ )]}' ' (H9) 

PHYSICAL LIBRATION, AND THE PRECESSION OF THE LUNAR EQUATOR 

Physical libration causes a periodic oscillation in the position of the lunar surface, thereby 
creating a small angular displacement between the principal lunar meridian and the earth-moon 
line of centers. The largest contribution to this displacement is given by 

X = a sin-Cg + (see Reference 4, p. 316) (120) 


where 


a = - 59" = - 2. 86 X 10’'' rad , 

Xo = constant dependent upon the Initial time, 


and 


= mean anomaly of the sun. 


Let c C and xy z be selenocentric, equatorial co- 
ordinate systems with the <f-axis directed toward the earth 
and the x-axis passing through the principal meridian (see 
Figure 7). Then, 


^ (?<f ’ ^ drj ’ 


X - (f cos X ^ V sin X > 


and 


/ 

© 


f = z 



X 


Figure 7 — Selenocentric, equatorial 
coordinate systems. 


-^siny + T 7 cos y . 


( 121 ) 
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Thus, 

X = cos X + sin X + 2arig cos-fig (tjcos X- ^ sin x) 

-ong^ sin-tg(T7Cosx-^ sinx) “ cos^ -fig (17 sin x + cos x) . (122) 

•where Og = mean motion of the sun's mean anomaly. But 

•i) cos X ~ t sin x “ '^Ug x cos -fig y 


and 


(9F 

^ cos X + •^ sin X - Sx 


(123) 


so Equation 122 becomes 


X = 2ang cos -fig (o-Og x cos -tg y) ~ sin -fig (i] cos x ‘S sin x) rig^ x cos 


.2 jf 


(124) 


If one neglects terms higher than third-order, Equation 124 reduces to 


^ + 2angy cos-fig . 


(125) 


A similar computation yields 


-9F . p 

- 2ong X cos -fig 


(126) 


Thus, 


in 

dx 


in 

(9y 


(127) 


where 


F = F 2aOg (xy - xy) cos -fig 


- F + 2an H cos -fig . 


(128) 
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Hence, the addition to the Hamiltonian of the problem is 


Then, 


and therefore 


FjClib.) - 2arig H cos . 

Fg^Clib.) - 0 

Faph'Clib.) - 0 

^3pr (lib.) = 2angHcosf.g , 


aSj' (lib. ) 

dT “ 2arigHc6s'C^. ^S^'Clib. ) = 2aHsin-t^. 


(129) 


(130) 


(131) 


Now, the precession in the lunar equator is created by the small angle of inclination ('d.'’32.'l) 
between the lunar equatorial plane and the ecliptic. This results in the regression of the equatorial 
node in the ecliptic. However, the motion of the node is already implicit in the hence no ad- 
ditional corrections need be made. 

COMPLETE "SECULAR" THIRD-ORDER HAMILTONIAH 

The complete secular third-order Hamiltonian is given by 

^3 ~ Fj" (coupling) + F3'(sun) + Fj" ( e) , (132) 

where F/ (coupling) is given by Equation 86 , F" (sun) by Equation 96, and F/ (a>) by Equation 117. 


THE DETERMINING FUNCTION FOR LONG-PERIOD TERMS 

Additional long-period perturbations depending on the motion of the node, the motion of the 
earth, and the motion of the sun will be given through the determining functions S^' and S 3 ' . The 
determining functions S 2 ' and S 3 ' are 

Sj' = S 2 ' (coupling) + Sj' (earth) (133) 
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and 


S3' = S3' (radiation) + S3' (sun) + S3' (e^) + S3' (ig) + S3' (®) + S3' (libration) ; (134) 

the various functions are defined by Equations 85, 92, 97, 102, 106, 111, 119, and 131. 

In the next section, the partial derivatives needed in order to find additional perturbations in 
the canonical elements are given. 


LONG-PERIOD PERTURBATIONS OF SECOND ORDER 


The perturbations of long period are obtained from 


dS 

as/ 

7 < 0. 

as/ 

as/ ^ 

T t 

as 

dh" 

L + 

aL" 

^ aL" 

aL" 

L 

■ az' ■ ^ 


as 

as/ 

as/ 



- ‘Ill , m 

ac" 

= g' + 

ac" 

^ ac" 

ac" 

G' 

ag' ^ ag' ag' ag' 


as 

as/ 

as/ 

as/ 

H' 

as 

as/ 

as/ 

^ <5S3' 

an" 

= h ' + 

aH" 

^ aH" 

^ aH" 

ah' 


ah' 

ah' 


(135) 


The terms corresponding to S/ have already been obtained. Next the partial derivatives of 
Sj' and S3' with respect to a", e", I", g', andh' are computed. For convenience, however, the 
primes have been dropped in this section. 


The "mean" mean motion n' (called n) depends on a' (called a) through the relation 





and 


as 

^ as 1 - as 


aL ^ 



as 

1 as cot I 

as 

ac ■ 

e y na de Jua {1 - e^] 

ai 
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and 


^ _ 1 ^ 

' sinl/Ma (l - e^) <51 

The following definitions are made: 

Cj = (l - e^) (2 “ 17e^) sin^ I cos I 

Cj = (2 + 3e^) ^ sin* I cos I 

C 3 = e*(l-e*](l + cosI)*(2“3cosI) 

= e*(l-e*)(l-cosI)*C2 + 3 cos I ) 

C 5 = e*sin*I(1 + cosI) [ 6 e*-(l + 5cosI)] 

Cj = e* sin* I( 1 - cos I ) [ 6 e * - ( 1 - 5 cos I )] 

Cj ' e*(l + cosI)* [2e * - ( 1 + cos I )] 

Cg = e* ( 1 - cos I)* [2e* - ( 1 - cos I)] 

Cg - (2 + 3e*) sin* I cos I 

Cjg = e*(l + cosI)* (l-'2cosI“5cos*l) 

Cjj - e*(l“COsI)* (1 - 2 cos I - 5 cog* l) 

Cj 2 = (2 + 3e *) COS I 

Cj 3 = e*(l + cosI)(3-5cosI) 

Cj^ = e*(l-cosI)(3+5cosI) 

Cj 5 = e* ( 1 + cos I ) ( 1 + 5 cos I ) 

Cjg = e* ( 1 - cos I ) ( 1 - 5 cos I ) 

Cj 2 = (1 + cos I ) [25 (1 - cos I - cos* 1-7 cos® l) - (1 - e*) (9 + 15 cos 1+15 cos* 1-175 cos® l)J 

Cjg' = (1-cosI) |^25(l+cosI-cos*I+7 cos® l) - (l-e*) (9-15cosI +15cos*I + 175cos®l)J 


(136) 
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Cjg = e^(l+cosI)^ (l + 2 cos I - 7 cos^ l) 

Cjo e^(l-cosl)^ (l - 2 cos I - 7 cos^ l) 

Cjj = (1+cosI) (3 - 25 cos I + 5 cos^ I + 25 cos^ l) 

C 22 “ (1-cosI) (3 + 25 cos I + 5 cos^ I - 25 cos^ l) 

C 23 = sin2 Icos I [ 35 - 10 ( 1 -e2) ( 2 + 7 cos=*l) - 3 ( 1 -6^)2 ( 1 - 14cos 2 l)J 

= e2 ( 1 + cos I)^ [35 (3 + 4 cos I - 30 cos2 I - 20cos^ I + 51 cos'* l) 

- 3 (1 - e^) (11 + 36 cos I - 126 cos^ I - 196 cos^ I + 315 cos'* l)j 

Cjj = e^(l-cosl)^ ^35 (3 - 4 cos I - 30 cos^ I + 20 cos^ I + 51 cos'* l) 

- 3 (1 - e^) (11 - 36 cos I - 126 cos^ I + 196 cos^ I + 315 cos'* l)j 

Cjg = e** sin^ I( 1 + cos I)^ (l + 2 cos 1-9 cos^ l) 

Cjj = sin^ 1(1 - cos I)^ (l - 2 cos I - 9 cos^ l) 

Cjg = (2 + 3e^) cos I (3 - 7 cos^ l) 

= e^(l + cos I) ( 3-19 cos 1-7 cos^ I + 35 cos ^ l) 

Cgg = (1 - cos I) (3 + 19 cos I - 7 cos^ I - 35 cos^ l) 

C 3 j = ( 1 + cos I) [35 (15 + 41 cos I - 262 cos^ I - 346 cos^ I + 519 cos'* I + 321 cos® l) 

+ 14 ( 1 - e"') (259 - 255 COS I - 566 cos^ I + 1710 cos® I - 525 cos^ I - 1455 cos® l) 

- (1 - e®)® (1559 - 1015 cos I - 5446 cos® I +7350 cos® I + 735 cos'* I - 5775 cos® l)j 

C 32 = (1 - cos I) [ 35 ( 15 -41 cos I - 262 cos® I + 346 cos® I + 519 cos'* I - 321 cos® l) 

+ 14 (1 - e®) (259 + 255 cos I - 566 cos® 1-1710 cos® I - 525 cos'* I + 1455 cos® l) 

- (1 - e®) ® (1559 + 1015 cos I - 5446 cos® I -7350 cos® I + 735 cos'* 1 + 5775 cos® l)] 
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( 1 + cos I)^ js (9 - 7 cos I “ 119 cos^ 1+15 cos® I + 174 cos'* l) 

- (1 - e®) (179 - 14 cos I - 1320 cos^ 1+30 cos® I + 1605 cos'* l) 

+ 3(l-e®)^ (I 8 + 7 cosI“135cos®I-15cos®I +165 cos'*!) 

(l~cosI)® j5(9+7cosI-119cos^I“15 cos® I + 174 cos'*!) 

- (1 - e®) (179 + 14 cos ! - 1320 cos® 1-30 cos® ! + 1605 cos'* l) 








+ 3(1 

-e®)® (18-1 

f cos ! - 135 cos® ! + 15 cos® ! + 165 cos'* 

Css 

= e® 

sin^ 

1(1 ■* 

cos!)® Jl 

- 9 cos® ! - (1 

-e®) (3-11 

cos® l)j 

C 36 

= e® 

sin ^ 

1(1 - 

■ CO s ! ) ® [^1 

- 9 CO s ® ! - (] 

1 - e®) ( 3 - 11 

cos® l)J 

^37 

-- (4 

+ 3e^ 

In 

+ cos ! ) ( 3 

- 5 cos 

!) (29 

- 126 cos® I ^ 

105 cos'* l) 

C 38 

= (4 

^ 3e^ 

0(1- 

-cos!) (3 

+ 5 cos 

I) (29 

- 126 cos® ! + 

105 cos'* !) 

C 3 9 

= e® 

sin^ 

1(1 ^ 

■ cos ! ) ( 3 ' 

- 5 cos 

I) ( 7 - 

15 cos® l) 


C 40 

- e® 

sin"^ 

1(1 - 

• cos I) (3 

+ 5 cos 

!) ( 7 - 

15 cos ® l) 



Then, the partial derivatives are as follows: 


rlSj' ( coupling) 
f9 a 



81926 


U'2 [32CjSi 


sin 2h - 4C, sin 4h + 80C, sin ( 2g + 2h) 


(1-e®)’ 

80C^ sin (2g - 2h) + IOC 5 sin (2g + 4h) + lOC^ sin(2g-4h) + 25Cy sin (4g + 4h) + 25Cg sin (4g- 4h) 


27 

* 

"‘e 


512e I „ * I J 2 


>(l-e®)* 


j^4Cg sin 2 h + 5Cj ^ sin ( 2 g + 2 h) + 5Cj ^ sin ( 2 g - 2 h)j 


27 2 n sin® ! 

10246 1 * J22^ 

/ 


,(l-e®) = 


8 Cj 2 (2 sin 2h - sin 4h) - 20Cj 3 sin (2g + 2h) 


( 137 ) 


-20Cj^ sin(2g-2h) +5Cjs sin(2g + 4h) +5Cjg sin(2g-4h)J 
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J2J2,b^ 


n sin^ I cos I 


(1-e^) 


117/n\2_,,^n sin^ I cos I . 

1 K't Sin 4 h 


7/2 


;in 2h + 32 ( *) ^22^ 




7/2 


Jab' ^ 2 ”^V-e^V ^ +C,g cos (g - 2 h) - 25 C ^9 cos ( 3 g + 2 h) - 250^0 cos ( 3 g - 2 h)] 


135 / n \ ^ j ne sin I 


315 

16384c 





(1-e^)^ 


8C23 sin 2h + 


sin ( 2g + 2h) + C,5 sin ( 2g - 2h) 


- 35C2J sin ( 4 g + 2 h) - 35C2, 


(4g-2h)] 




n sin^ I 
( 1 - 02 ) 11/2 


ne s in I _ r 
a"(l-02)5 L^31 


[^C,s 


sin 2h + Cjg 


cos (g + 2 h) + C32 


sin (2g + 2h) + Cjp sin ( 2g - 2h)j 


cos ( g - 2h) - I4C33 cos ( 3 g + 2h) 


- I4C34 cos ( 3 g - 2 h) + IO5C35 cos ( 5 g + 2 h) + IO5C35 


cos 


(5g-2h)] 


855 
" 8192 



ne si n I 

a^(l-e2)l'’^ 


1^X3, cos (g + 2h) + 2033 cos (g - 2h) 


-yCjg cos (3g + 2h) -7C^o cos (3g-2h)] • (138) 


ciSj' (coupling) 

s'l 



na^ sin I 


16 (1 - 3 cos2 l) (1 - e^) (2 - 17e2) sin 2h 


+ 2 (1 - 3 cos 2 l) (2 / 3e2 )2 sin 4h - 40( 1 + cos I ) ( 1 - 9 cos I ) (l - e^) sin ( 2g + 2h) 

+ 40( 1 - cos I) ( 1 9 cos I) 02 (1 - 02) sin (2g - 2h) 10( 1 + cos I) e2 1^/3 - 5 cos I) ( 1 + 2 cos I) 


- 3e2 ( 1 - 3 cos I) sin (2g + 4h) - 10( 1 - cos I)e2 |^(3 + 5 cos I) ( 1 - 2 cos I ) 
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- 3e ^(1 + 3 CO si)] sin(2g-4h) + 25(1 +c'os I)*e^ [2( 1 + cos I ) - 3e^] sin ( 4g + 4h) 


25( l-cosI)^e^ [2(l-cosI)-3e^] sin(4g-4h)^ 


9 / ”e 


^ 2 ^^ 2 ^ [' (l-3cos^l] (2+3e^) sin 2h 


5( 1 + cos I ) (l - cos I - 5 cos^ l] sin (2g + 2h) + 5( 1 ~ cos I) (l + cos 1-5 cos^ l] sin (2g - 2h)^ 


/ n \ 

9 / e\ . ns in If, ■, \ / „\ 

256e ( *) J 22 '^^ /i _ 2\2 L4(1-3cos^I) (2 + 3e^) ( 2 sin 2h - sin 4h) 

\"e/ r ^ ) 


20( 1 + cos I) + 7 cos I - 10 cos^ l) sin (2g + 2h) 


20( 1 - cos I ) (l ■ 7 cos I - 10 cos^ l) sin ( 2g - 2h) 


5(1 + cos I) (1 + 2 cos I) (3-5 cos I) sin(2g + 4h) 


5( 1 - cos I ) ( 1 - 2 cos I ) ( 3 + 5 cos I ) s in ( 2g - 4h)J 


9 ^_n^2 
In 


E , 


n sin I (l - 3 cos^ l) 
(l - e' 


8 I *) J2J22*^'' ‘ sin2h + 


9 

i n„ 


16l„*l J2'2b" 




n sin I (l - 3 cos^ l) 
a2 .^2)77^ 


sin 4h 


1024s I , * / J 3 
"e 




(l-e = )' 


{(1+cosI) [25cosl(5 + 19 cos I 1 3 cos^ 1-35 cos^ l) 


+ (1 - ) (24 + 27 cos I - 555 cos ^ I - 235 cos ^ I + 875 cos'’ l)J cos ( g + 2h) 

+ ( 1 - cos I ) |25 cos I (5 - 19 cos 1 + 3 cos ^1 + 35 cos^ l) 

- (1 - e^) (24 - 27 cos I - 555 cos^ I + 235 cos^ 1 +875 cos'’ I ]j cos (g - 2h) 

+ 25e^(l+cosI)^ (4-13cosI-22cos^I+35cos®lj cos(3g + 2h) 

- 25e^ ( 1 “ cos I ( 4+13 cos 1-22 cos^ 1-35 cos^ l] cos (3g - 2h)^ 
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9 / n 5 ne 

128 J 3 J22^ 


j( 1 + cos I ) (22 + 21cosI-155cos^I~5cos^I + 125cos'’l) cos(g+2h) 


45 4 n sin I 

J 4 ^ 


40966 \^n;y 


( 1 -e^)' 


-(1-cosI) ^22 - 21 cos I - 155 cos ^ I + 5cos^ I + 125 cos'* l) cos ( g ~ 2h)j 
-^4 jsS (1 - 3 cos^ l) - 10 (1 - e^) (2 + 15 cos^ I ~ 35 cos'* l) 


-3(1- 

e®) ^ ( 1-45 cos® I 

-3(1- 

e®) (29 - 72 cos I - 

- e® (1 

- cos I) [ 35(5 + 24 

-3(1- 

e®) (29 + 72 cos I - 

- 35e‘* 

[(l + cos I)® (2-9< 


- ( 1 - cos I (2 + 9 cos 1-14 cos^ 1 + 27 cos^ l) sin (4g - 2h) 


JU (t^ ) J 4 J 22 7 T|f- 7 ^ilT 2 - [4 (2 ^ 3e^) (3 - 30 cos^ I + 35 cos^ l) sin 2 h 


2e^ ( 1 + cos I) (8 + 21cosI-87cos^I-35cos^I + 105cos‘*l) sin(2g+2h) 

+ 2e^ ( 1 “ cos I) (8 - 21 cos 1-87 cos^ 1+35 cos^ I + 105 cos'* l) sin ( 2g - 2h)J 


45 / **e \ 5 n e_ _ 

+ 327686 Jsb ^3(i_;2-)5 

- 2793 cos® I - 2247 cos® l) + 14 (l - e^) (4 - 1905 cos I + 5329 cos^ I + 1874 cos® I 

- 16350 cos'* I + 1695 cos® I + 10185 cos® l) - (l - (544 - 15025 cos I + 19649 cos® I 

+ 32074 cos® I - 64890 cos'* I - 10185 cos ® I + 40425 cos ® l)J cos ( g + 2h) 


{■ 


-(1 + cosI) 


[^35 (56 - 


513 cos I - 1423 cos® I + 2778 cos® I + 3854 cos'* I 
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+ (1 - cos I) [ 35(56 + 513 cos I - 1423 cos^ I - 2778 cos ^ I 

+ 3854 cos'* I + 2793 cos® I - 2247 cos® l)+ 14(l - e*) (4 + 1905 cos I + 5329 cos^ I - 1874 cos® I 

- 16350 cos'* I - 1695 cos® I + 10185 cos® l) - (l - (544 + 15025 cos I + 19649 cos® I -32074 cos® I 

- 64890 cos'* I + 10185 cos® I + 40425 cos® cos (g - 2h) 

+ 14(1 + cos I)® [5 (11 - 279 cos I - 165 cos® I + 1321 cos® I + 258 cos'* I - 1218 cos® l) 

- (1 - e®) (344 - 3205 cos I - 2494 cos® I + 13080 cos® I + 3030 cos'* I - 11235 cos® l) 

+ 3 ( 1 - 6 ®)® (43 - 310 cos I - 343 cos® I +1305 cos®I + 420 cos'* I - 1155 cos® l)J cos ( 3g + 2h) 

- 14(1 - cos I)® [5 (11 + 279 cos I - 165 cos® I - 1321 cos® I + 258 cos'* I + 1218 cos® l) 

- (1 - e®) (344 +3205 cos I - 2494 cos® I - 13080 cos® I +3030 cos'* I + 11235 cos® l) 

+ 3(1-0®)® (43 + 310cosl - 343 cos®I - 1305 cos®I + 420 cos'* 1 + 1155 cos®l)j cos(3g~2h) 

- 105e® sin® 1(1+ cosi)® [ 2 - 23 cos I - 18 cos® I + 63 cos® I 

- (1 - e®) ( 6-37 cos 1-22 cos® 1+77 cos® l)J cos ( 5g + 2h) 

+ 105e® sin® 1(1- cos I)® [ 2 + 23 cos I - 18 cos® I -63 cos® I 

- (1 - e®) (6 + 37 cos 0-22 cos® 1-77 cos® l)j cos (5g - 2h) 

+ 4096 ^;;^^ "J^y777^|2(4 + 3e®)[(l+cosI) (58 + 1075 cosI- 1947 cos2l 

- 3402 cos® I + 5460 cos'* I + 2415 cos® I - 3675 cos® l) cos (g+ 2h) 
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(1 - cos I) (58- 1075 cos I - 1947 cos^ I + 3402 cos^ I 


+ 5460 cos'* I - 2415 cos® I - 3675 cos® l) cos ( g - 2h)J 
- 7e^ sin^ I [^(1 + cos I) (l4 + 209 cos I - 355 cos^ I 


345 cos® I + 525 cos'* l) cos ( 3g + 2h) - ( 1 - cos I) (l4 - 209 cos I - 355 cos® I 


dS' ( coupling) 


na®e 


d e 


4096^2 



(l-e®)^ 


345 cos^ 1 + 525 cos'^ l) cos ( 3g - 2h) 1 I , (139) 


96 sin® I cos l(l2-17e®) (l-e®)sin2h 


- 4 sin® I cos I (2 + 3e®) (l4 -9e®) sin 4h + 80( 1 + cos I)® (2 - 3 cos I) (2 - 3e®) (l - e®) sin (2g + 2h) 
+ 80( 1 - cos I)® (2 + 3cos I) (2 - 3e®) (l - e®) sin (2g- 2h) 

+ 10 sin® 1(1 + cos I) [oe® (4 - 3e®) - (2 - e®) ( 1 + 5 cos I)j sin ( 2g + 4h) 

+ 10 sin® I( 1 - cos I) [^6e® (4 - 3e®) - (2 - e®) (1-5 cos I)j sin (2g - 4h) 

+ 25( 1 + cos I)® j^2e® (4 - 3e®) - (2 - e®} ( 1 + cos I )j sin (4g + 4h) 

+ 25( 1 - cos I)® [^2e® (4 - 3e ®) - (2 - e®) ( 1 - cos I)j sin ( 4g - 4h)|- 


9 / "(E \ , , ne 

+ 1 no ^ ( * 1 J ^ b 


J2 


(l-e®)^ 


1^4 sin® I cos 1 ( 7 + 3e®) sin 2h 


5(l + cosI)® (l + 2cosl-5cos®l)(l + e®) sin(2g + 2h) 


5(l-cosI)® (l-2cosI-5cos®lj (l+e®)sin(2g-2h) 


68 



y / t 2 s irv^ I 

It 0 \ ^ 

\"(D / 




{ 1 -e^Y 


j^- 8 cos I (7 + 3e^) (2 sin 2h - sin 4h) 


-20(l + cosl)(3-5cosl) (l+e^) sin(2g+2h) 


- 20( 1 - cos I) (3 + 5 cos I ) (1 + e^) sin (2g- 2h) 


5( 1 + cos I) ( 1 + 5 cos I) (1 + e^) sin ( 2g + 4h) 


5(1“ CO si) (l“ 5 cos I) (1 + e^) sin(2g~4h)J 


63 / n 


ne sin^ I cos I 


/ \ O / **\ o A 


63 / n\^ „ ^ ne sin^ I cos I 

I 1 T h ’ 

16 ^n^y J 22 g2^j^_p2^9/2 sin 4h 


I- r^y ,3 

1024 e(„*l J3'’ 


n sin I 


f( 1 * cos I ) ^150 (^1 - cos I - cos^ I - 7 cos^ l) 


“ (l61-65cosl-65cos^l“1575 cos^l) + 3(l-e^)^ ( 9 * IScosI 

+ 15 cos^ I - 175 cosM)] cos (g H 2h) + ( 1 - cos I) [l50(l <■ cos I - cos ^ I <■ 7'cos^ l) 
■ (l“e^) (I 6 I '65cosI“65cos^I< 1575 cos^l) 

• 3(l“e^)^ (9“15cosI ’ 15cos^I • 175 cos^ cos ( g - 2h ) 

- 75e^ (1 » e^) j^( 1 * cos I ) ^ ^1 * 2 cos I - 7 cos^ l) cos ( 3g + 2h) 


( 1 - cos I ) ^ 2 I “ 7 cos ^ l) cos ( 3g - 2h) 


9 / n \2 sin I 

^ 128 j J3 J22 (1 - [( 1 + cos I ) (3 “ 25 cos I + 5 cos^ I 


25 cos ^ l) cos (gT2h) + (l“cosI) (3+25 cos 1 + 5 cos^ I - 25 cos ^ l) cos ( g - 2h)] 


19 


il- 


+ {ssin^lcosj [70- 15(1.- e^) (2 + 7cos2l) 

-3(1-6^]^ (l-14cos^l)J sin 2h 

+ ( 1 + cos I ) ^ J^70 (3+4 cos 1-30 cos ^ I - 20 cos^ 1 + 51 cos'* l) 

- 3 (1 - e^) (69 + 124 CO si - 714cos^I - 644 cos^ I + 1365 cos'* l) 

+ 3(l~e^)^ (ll + 36cosI-126cos^I-196cos^I+315cos'*l)J sin(2g + 2h) 
+ (1 - cos 1)2 [^70 (3 - 4 cos I - 30 cos^ 1 + 20 cos^ 1 + 51 cos'* l) 

- 3 (1 - e2) (69 - 124 cos I - 714 cos^ I + 644 cos^ I + 1365 cos'* l) 


+ 3 (l-e 2)2 ( 11-35 cos I - 126 cos^ I + 196 cos^ I + 315 cos'* l)j sin ( 2g - 2h) 
- 35 e 2 (i+e2)sin2l ^(l+cosI)2 (l + 2cosI-9cos2l) sin(4g+2h) 


(1- cos 1)2 (1 - 2 cos I - 9 cos2 I^ sin(4g-2hij^ 


45 / n' 


J 4 J 22 >^® {'^(28 " 27e2) cosl(3-7 cos2l) sir 


+ (2 + 9e2) (1 + cosI) (3-19cosI-7cos2l + 35cos^l)sin(2g-'2h) 


(1-cosI) (3+19cosI“7cos2l-35cos^l) sin(2g- 


45 / \ c n sin I ( r , „ 

+ 327^6 ( * ) J5 3 , jw i( 1 + cos I) 1^350 (l5 + 41 cos I - 262 cos2 I 

V+'e / ® ^ / I 

- 346 cos^ I + 519 cos'* I + 321 cos^ l) 

+ 7 (1 - e2) (3469 - 5925 cos I + 2734 cos2 I 
■ + 42930 cos^ I - 31755 cos'* I - 37725 cos® l) 


70 



8 (l - 0^)2 (4342 - 3885 cos I - 11018 cos^ I + 26460 cos^ I 


5880 cos'* I - 22155 cos= l) + 5(l - (l559 - 1015 cos I - 5446 cos^ I 


7350 cos^ I + 735 cos'* I - 5775 cos = 


»] 


cos ( g + 2h) 


( 1 - CO s I ) 


[ 350(15 - 41 COS I - 262 cos^ I + 346 cos^ I + 519 cos'* I - 321 cos= l) 


7 (1 - 0 = ) (3469 + 5925 cos I + 2734 cos^ I - 42930 cos^ I 


31755 cos'^ I + 37725 cos= l) - 8(l - 0^)2 (4342 + 3885 cos I 


11018 cos^ I - 26460 cos^ I - 5880 cos'* I ♦ 22155 cos^ l) 


5 (1 - Y (1559 * 1015 cos I - 5446 cos^ I - 7350 cos^ I 


735 cos'* I f 5775 CO s^ l] cos (g~ 2h) 


14( 1 * cos I )* 


50 ( 9-7 cos 1-119 cos*' I < 15 cos*' I » 174 cos'* 


(1 - (1837 - 427 cos I - 15915 cos^ I ■ 915 cos' I • 20670 cos'* l) 


(1 - e')' (1577 + 28 cos I - 11670 cos' 1-60 cos' I + 14205 cos'* l) 


15 ( 1 - 0 ']' (l8+7cosl-135cos'l-15cos'l + 165cos'*l)j cos(3g 


2h) 


14(1 - cos I)* 


[50 (9 + 7 cos I - 119 cos' I - 15 cos' 1 + 174 cos I ) 


(1 - e') (1837 + 427 cos I - 15915 cos' I - 915 cos' I + 20670 cos'* l) 


(1 - e')' (1577 - 28 cos I - 11670 cos' 1 + 60 cos' I + 14205 cos'* l) 


71 


- lS(l“e^)^ (l8 - 7 cos I " 135 CO I + 15cos^I + 165 cos^ l)j cos ( 3g - 2h) 

+ 105e^ sin ^ I [ 10 ( 1-9 cos^ l) - (1 -e2) ( 31 - ISlcos^ l) 

+ 5 ( 1 - (3 - 11 cos^ ■)][(> + cos I)^ cos ( 5g + 2h) + ( 1 “ cos I)^ cos ( 5g - 2h)J^ 


+ Jg ‘TTff 'jUs" ; -[2(4 + 5762 +30e'*) (29 - 126 cos2l 

C / ^ ® ) L 


4096 


. n 


+ 105 cos ^ l) [( ^ ^ I ) ( 3 - 5 cos I ) cos ( g + 2h) + ( 1 - cos I) (3 + 5cosI) cos ( g - 2h )J 
-7e2 (3+10e2) sin2l(7-15cos2lj [('l + cosI)(3-5 cos I) cos (3g + 2h) 

+ ( 1 - cos I ) ( 3 + 5 cos I ) cos ( 3g - 2h")j^ , (140) 


dSj' (coupling) 
f5g 



(l-e2) 


^ ^ [8C3 cos (2g + 2h) + 8C4 cos (2g - 2h) 


■ Cj cos ( 2g + 4h) + Cg cos ( 2g - 4h ) ^ 5Cj cos ( 4g + 4h ) + 5Cg cos ( 4g - 4h)J 


45 

128e 


J,b2 


(l-e2)^ 


[Cjg cos C2g + 2h) + Cj j cos ( 2g - 2h)J 



n sin^ I 

(T^)2 


[-4Ci3 cos (2g + 2h) -4 Cj4 


cos ( 2g - 2h) 


+ Cj 5 cos (2g +4h) + Cjg cos(2g-4h)J 


/n \ 2 

9 / E \ ne sin I 

a(l-e2)'3 


[Cj, sin (g + 2h) + Cjg sin ( g - 2h) 


“75Cjg sin(3g-f 2h) - 75 C 2 Q sin(3g-2h) 


9 « / n \ 2 ^ ne^i n JI 

128 \^*j J3J22b ^3 ^i_g'2)9/2 


C 21 sin (g + 2 h) + C 22 sin (g - 2 h)j 


72 



45 

4096 e 





(l-e=V 


cos ( 2 g + 2 h) +C 25 cos ( 2 g - 2 h) 


- 70C2g (4g + 2h) - 7 OC 2 , cos (4g - 2h)] 


, '^5 / n ^ ^ , fi n sin^ I r -1 

256 j J 4 J 22 b g 4 ^i_g 2 yi /2 [C 29 cos( 2 g + 2 h) +C 3 t,cos( 2 g- 2 h)J 

/ \ 2 

45 [ \ _ ne sin I r 

32768f J 5 ^ ^3 [_^3 i ^in (g + 2h) +C 32 sin (g - 2h) - 42 C 33 sin (3g + 2h) 


- 42 C 34 sin (3g - 2h) + 525 C 35 sin (5g + 2h) + 525C3g (5g " 2h)j 

45 / n_\ ^ , ne s in I r 

4096 f*j JsJ 22 b g 5 ^i_p 2 ji 32 [ 2 C 37 sin (g 2h) t 2033 sin (g - 2h) 

- 2 IC 3 , sin (3g t 2h) - 2 IC 40 sin (3g- 2h)] , ( 141 ) 

^Sj’Ccoupl^ _ _9 na^ r 

(9h 1024f^y^y 2h 4 C 2 co s 4h + 4 OC 3 cos ( 2g ^ 2h ) 

- 40C^ cos ( 2g - 2h) + IOC 5 cos ( 2g 4h) - lOCg cos ( 2g - 4h) 


+ 25C, cos ( 4g + 4h) - 25Cg cos ( 4g - 4h 



4C„ cos 2h + 5C, n cos 


(1 - e^) 2 L '"-9 ""'lo 


( 2 g + 2 h) - 5Cj j cos ( 2 g - 2 h)J 




n sin ^ I 


[" 2 2h - cos 4h) - IOC , 3 cos ( 2g + 2h) + lOC,^ cos ( 2g - 2 h) 


+ 5 C ,5 cos (2g + 4h) - 5C,g cos ( 2g - 4h)j 


73 



sin^ I cos I 




ne sin I 
a(l-e2) = 


COS 2h - "7 


9 / n 


7/2 4 1 „ * J J22 „2 


n sin^ I cos I 


cos 4h 


[Cj7 sin(g + 2h) -Cja sin(g- 2h) - 25Cjg sin (3g + 2h) + sin (3g- 2h)] 


9 / n \ ^ 5 ne si n I 

64 J 3 J 226 g3^i_^2^9/2' 


jCjj sin (g + 2h) - sin ( g - 2h)j 


45 

4096 e 



J4b^ 


(1-e^)' 


8C23 COS 2h t C24 cos ( 2g + 2h) " C25 cos ( 2g - 2h) 


- SSCjg cos (4g - 2h) * 35C2, cos ( 4g - 2h)j 


45 ! n^^ 
256 

■ ‘e , 


n\* nsin^Ir H 

^ J 4 J 22 ^^ g4^^_p2^ii 2 j^4C2a cos 2h Cjg COS ( 2g + 2h) - C 3 Q cos ( 2g 2h )J 


45 / "e 


J5 [Cjj sin (g - 2h) - C32 sin (g - 2h) - I4C33 sin ( 3g " 2h') 


16384£ 


(1-e^) 


HCj^ sin (3g - 2h) + IO5C35 sin ( 5g + 2h) - 105Cjg sin ( 5g “ 21i)j 


45 / n 


[2C37 sin (g4 2h) -2C33 sin(g-2hi 


2048 


■ yCjg sin ( 3g + 2h) - 7C^g sin ( 3g - 2h)J , (142) 


riSj' ( earth) 
d a 


512e 


— % 1 r~ 9( 1 ■* cos I) (l-^ 10 CO si" 15 cos (4^3e^) sin(g+h) 

/ ^(£ L ^ ^ 


+ 9(1- cos I) (l-10cosl-15cos^l) (4+3e^)sin(g-h) 


+ 15sin^I(l ^cos I) (4 +3e^) s i n ( g + 3h ) 


-15sin^I(l-cosI) (4 + 3e^) sin(g“3h) 
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+ 315 sin^ 1(1 +COS I)e^ sin(3g+h) 


i^Sj' (earth) 
51 


dS^ (earth) 
de 


-315 sin^ 1(1 - cos I)e^ sin(3g-h) 

+ 35(1 +COS 1)3 sin(3g + 3h) -35(1 - cos I)® sin (3g - 3h)J , (143) 
sin I 

p (ll - 10 cos I - 45 cos^ l) (4 + 3e3) sin(g+h) 

+ 3 (11 + 10 CO si -45cos3i) (4 + 3e^) sin(g-h) 

- 5( 1 + cos I ) ( 1 - 3 cos I ) (4 + 3e3) sin ( g + 3h) 

- 5( 1 - cos I)(l+3cosI) (4 + 3e3) sin ( g - 3h) 

- 105(1 +COS I) (1 - 3 cos I) e^sin (3g + h) 

- 105( 1 - cos I ) ( 1 +3 cos I) sin (3g-h) 

- 35( 1 + cos I)3e3 sin(3g + 3h) - 35( l-cosI)3e3sin(3g-3h)j , (144) 

^-3(1 + cos I) (1 + 10 cos 1-15 cos^ l) (4 + Oe^) sin (g + h) 

+ 3(1- CO si) (1-10 cosI-15cos3i) (4 + 9e3)sin(g-h) 

+ 5sin3l(l + cosI) ^4 + 9e3) sin(g + 3h) 

- 5 sin^ I( I - cos I) (4+9e3) sin(g-3h) 

+ 315sin3l(l +COS I)e3 sin(3g+h) 

-315 sin^ 1(1 - cos 1)6^ sin(3g-h) 

+ 35( 1 +cosI)3e3sin(3g + 3h) - 35( 1 - cos I)3e3sin(3g-3h)J , (145) 


5 ^ 

512e ( * ) a 
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[- 3(1 + cos I) (l + 10 cos 1-15 cos^ l) (4 + 3e^) cos (g + h) 

+ 3(1 -cos I) (l - 10 cos I - 15 cos^ l) (4 + 3e^) cos ( g - h) 

+ 5 sin^ 1(1 + cos I) (4 + 3e^) cos(g + 3h) 

- 5 sin^ 1(1 - cos I) (4 + 3e ^)cos(g-3h) 

+ 315 sin ^1(1 +COS I) cos(3g+h) 

- 315 sin^ 1(1 - cos I) cos (3g-h) 

+ 35(1 + cos I)^ cos (3g + 3h) - 35( 1 - cos I)^ cos ( 3g 


dS^ (earth) 
5g 



5Sj' (earth) 

5h 


[■ 3 (l + cos l) (l + 10 cos 1-15 cos^ l) (4 + 3e^) cos (g + h) 

- 3( 1 - cos I) (1-10 cos 1-15 cos^ l) (4 + 3e^) cos (g- h) 

+ 15 sin^ I( 1 + cos I) (4 + 3e^) cos (g + 3h) 

+ 15 sin^ 1(1 -cos I) (4 + 3e^) cos (g -3h) 

+ 105 sin^ 1(1 +COS I) cos(3g + h) 

+ 105 sin^ 1(1- cos I) cos (3g -h) 



5 S 3 ' ( radiation) 
da 


(5 S 3 ' (radiation) 


+ 35( 1 + cos I)^ cos (3g + 3h) + 35( 1 - cos I)^ cos ( 3g 


3 ere 


9 * * 

"e “"e 


i sin(h + g + + sin^ ^ sin(h - g + - \ J 


3 crae sin I 

~4 * * 


sin (h + g + + sin (h - g + . 


3h)] , ( 146 ) 


3h)] . ( 147 ) 

( 148 ) 

( 149 ) 
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( radiation) 

de 


cos 2 2 sin(h + g + \^-Xj +sin 2 ^ sin (h - g + 
"e "ffi L 


i9Sj' (radiation) 

5g 


3 crae 


2 * * 
"e '^e 


cos^ ~2 cos (h + g + - sin^ “2 c°s (h 


g + 




< 983 ' ( radi ation) 2 


dh 


7 * * 


cos 2 2 cos(h + g+X^-Xj +sin 2 i cos (h - g + X^ - X^) 


i3 S 3 ' ( sun ) 2 ^ 

^ [2(2 + 3e^) sin^I sin2(h + X^-Xg) 


d a 


+ 5e^ ( 1 + cos I sin 2(h + g + X^ - X^^) + 5e^ ( 1 - cos I) ^ sin 2 (h - g + X^ 


(9S3'(sun) 3 n^^a^sinl 

51 " 32 „ * [ 2 ( 2 +3e2) cosi sin 2(h +X^ -Xg) 

I * « 1 I _ U 


“ 5e^ ( 1 + cos I ) sin 2 {h + g + X^ - X^) + 5e^ ( 1 “ cos I) sin 2(h - g + X^ 


5S,' (sun) 


3 ni e 


offi r ^ , 1 

5e " 32„*_„* 6 sin-^ I sin 2(h +X^-Xg) 


+ 5(1+ cosl)^ sin2(h + g+ X„ - X^) + 5( 1 - cos I)^sin2(h-g + Xg 


5 S 3 ' ( sun) 

5g 


15 
32 


* _ * T( 1 + cos I cos 2 (h + g + X^ - Xg) - ( 1 - cos I)^ cos 2 (h - g + Xgj - X 


dSj ( sun) nj 


5h 


32 n*°-n* [2 (2 + 3e^) sin^ I cos 2 (h + X^ - X^) 


5e^(l + cosI)^cos2(h + g+ X^-Xg) +5e^(l-cosI)^cos2(h-g+X^ 


(e^) 


S a 


32e "a; ■ 


j“ 2 (2 + 3e^) (1 “ 3 cos^ l) + 30e^ sin^ I cos 2g 


(150) 

(151) 

(152) 

-K)] ■ (153) 

-^o)] • (154) 

-K)] ’( 155 ) 
. 3 )] . (156) 

-^o)] .(157) 
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+ 6 (2 + 3e^) sin^ I cos 2h + 15e^ ( 1 cos I)^ cos (2g + 2h) 


^53' (ej 


ai 


'^^3 (^e) 

de 


n_ e„ a‘ e 


32e C 


‘’S3' (eg) 45 


<5g 


‘’S 3 ' K) 9 


dh 


+ 15e^ ( 1 - cos I) ^ cos ( 2g “ 2h)J sin I ^ 
3 ^ rijj Cg, sin I 2 (2 + 3e^) cos I + lOe^ cos I cos 2g 

2 (2 + 3e^) cos I cos 2h - 5e^ ( 1 + cos I ) cos ( 2g + 2h) 

+ 5e^ ( 1 ” cos I ) cos ( 2g - 2h)J sin , 

2( 1 “ 3 cos^ I ) + 10 sin^ I cos 2g + 6 sin^ I cos 2h 

S(l + cosl)^cos(2g + 2h) + 5(l-cosI)^cos(2g~2h)j sin 

‘4'd r 

3 ^ n^je^a^e^ 2 sin^Isin 2 g + (l+ cos I)^sin( 2 g+ 2 h) 

(l-cosI)^sin(2g-2h)j sin 
3 ^ n^Cija^ j2(2 + 3e^) sin^Isin2h+5e^(l+cosI)^sin(2g+2h) 

(l-cosI)^sin(2g~2h)J sin I ^ 


Se-^ 


^S3 (V) 


3 ng2 a sin I sini^ 

d a 


- 4e + [- 2 (2 + 3e ) 


‘’S3' (V ) 


+ 5e^(l+cosI)sin(2g+h) + 5e^(l“cosI)sin(2g“h)J 
[2(2 + 3e2)(l-2cos2l) sinh 


^ n^a^ sini 
8e 


“5e^(l +COS I)(l “2 CO si) sin(2g+h) +5e^ (l-cosI)(l + 2cosI) sin(2g~h)J cos 


(158) 


(159) 


(160) 


(161) 


(162) 


, (163) 


(164) 
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^^3 (^e) 

(9 e 


3 e sin I sin i^ 


■ 6 cos I sin h + 5( 1 + cos I)sin(2g + h) 


5( 1 - cos I) sin ( 2g - h)j cos , (165) 


^^3 _ 15 ^2 sin I sin ij 


[^( 1 + cos I) cos (2g + h) + ( 1 - cos I) cos ( 2g -h)J 


(ie) 


a2 sin I sini^^ 


2 (2 + 3e2) cos I cos h 


+ 5e2 ( 1 + cos I) cos (2g + h) - 5e2 ( 1 - cos I) cos ( 2g - h)J cos , (167) 


'^S' (*) 


rlSj' ^ ^ H(j 2 j 2 R^2 

f5 1 16 e „ * _ nj 

"c '^Cc 


•|- 2 l^sin I cos I sin sin 2(h + - Cl') 


4 (1 - 2 cos2 l) cos sin (h + - fi')j 

(l + 2yT-e2) ^sinI(l + cosI)sinI,j^sin2(h+g+X^-r;' 


sinl(l - cos I) sin I,j,sin 2(h~g + Q') 


4(l + cosI)(l-2 cos I ) cos Igj sin (h + 2g + - 0.') 


4(1 - cos I) (1 + 2cos I) cos Igj sin (h - 2g + Xgj- O')]} . (169) 


ciSj' (e) 
d e 


j 2 K,? /3^ sin (2 + yf-e2) 
e(n;-NnJ 


j^- ( 1 + cos 1)2 sin sin 2(h + g + - O') 


(l-cosI)2 sinl^ sin2(h-g+\^g-0') 
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(e) 

i3g 


aSj' (®) 

ah 


+ 8 sin 1(1 + cos I ) cos sin (h + 2g + - Cl') 


8 sin I( 1 - cos I) cos sin (h - 2g + - n')J , (170) 


( 1 + cos I sin cos 2 (h + g 


16e 


^Oc) 


+ \ -n' 


( 1 - cos I ) ^ sin Ig cos 2 (h - g + - Cl') 

8 sin 1(1 + cos I) cos cos (h + 2g + ~ O') 


+ 8 sin I( 1 - cos I ) cos Ig cos [h - 2g + - n')J , (171) 


„ 2 ; D 2 
16« * 


3 




2 sin I jsin I sin cos 2(h + ^^ “ Cl') 


+ 4 cos I cos I . cos 


(h + X.^ - n')] 

/5^ (l + 2 Vl - e^) |^-(l + cosI)^ sin cos 2(h + g + X^ - Cl') 


- ( 1 - cos I ) ^ sin cos 2 (h - g + - f!') 


+ 4 sin 1(1+ cos I ) cos cos (h + 2g + X^ - Cl') 


4 sin I( 1 ~ cos I) cos cos (h - 2g + X^ - ’ (1'^^) 


aS3 (lib. ) 

da 


0 , 


(173) 


as/ (lib. ) 
ai 


2a G sin I sin 


(174) 


as/ (lib. ) 

de 


= 0 


(175) 


as/ (lib. ) 
ag 


0 . 


(176) 
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and 


dS' (lib. ) 

ah "" 0 • 


(177) 


SECULAR PERTURBATIONS AND PERTURBATIONS DEPENDING STRICTLY ON g" 

At this stage of the problem, the remaining part of the Hamiltonian is 


3 


3; . 3; 


3 '' 


. 3; 


(178) 


where 


and 


3o" 


const . 


2L'2 

3^" - H ' = const . , 

^ F;(-, b", L', G", H“) 


3," -- f;(-, g". L', G". H ) 


andL' , T, and H" are constants in time. T is effectively a third-order quantity; but in the manipu- 
lation of von Zeipel's method its contributions are of the zero order, for it is a variable itself, 
independent of all the others. 

It is easy to see that the variable part of the Hamiltonian which will generate the equations 
of motion, is factored by small parameters. Therefore, a method of successive approximations 
such as von Zeipel's cannot be applied. The system of differential equations produced by 3" must 
be integrated directly or by using some other kind of approximation. 

In a previous version of this theory (Reference 5), the system of differential equations 


G" 


.IF" 

<5g" 


1“ 


<9L' 


dF" 

.9G" 
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and 


h" = 


ar" _ 
an" 


where F" = Sq" + F^' +F 2 ", was integrated by making use of a method involving elliptic integrals. 
This method has appeared in Reference 6, where it was applied to a different problem. In the 
previous version, it was assumed that J3, J4, J5, and all other spherical harmonics of the moon 
were of third order or smaller, making the method involving elliptic integrals applicable. Since 
these and other spherical harmonics of the moon are probably not that small (Reference 3, for 
example), this method is no longer applicable. Thus it is suggested that the equations of motion 
for the secular Hamiltonian F" be integrated numerically. We have 


F" = 


p " + F " + F " + F " 

^0 2 3 


where 


Fo" 


2L' 


f; 


He H" 


(179) 


Fj" is given by Equation 78, and Fj" is given by Equation 132. The system of equations to be inte- 
grated numerically is as follows: 


^ dji ^ 

i9g" (9g" 


.9F" 

(9L" 


aV <9L" 


^F" 

_ ajl 

clF 

dC 

<?G" 

<9G 


and 


.9F" 

♦ 2 

aF 

aw 

“ n- " 

® SU" 

aw 


(180) 
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Again, it is more convenient to compute the partials with respect to Keplerian elements, and use 
the expressions 


aF" 

aL" 


2 /il ^ ^ ^ aF" 

V ^ a a ” e " Vua” a e “ 


aF" 

ac" 


J_ /l - e" 
e" y /za" 


aF" ^ cot aF^ 

(1 ai" 


and 


aF" 

an" 


ap" 


sin I " y^a ■ (l - e " 


ii^ai' 


Thus, 


aPj" apj" 1 ( VJj- r 

^ ' aT^ " ^ (2* (1 " 3 cos2 I") 4 15e"2 sin^ I" cos 2g'] 


1u2t M (l - 3 cos^ I") 3 A<e" sin I" (l - 5 cos^ I") 


+ 4 t>^J 


a"4 (i -e"2)3. 


2 


a'-(l-e-)^ 


sin g 


a"«(l-e"2j^ [(2 + 3e"2) (3 - 30 cos ^ I " + 35 cos '' I ") 


- lOe"^ sin 2 


I" (l - 7 cos^ I") cos 2g ] " b5 Jj [b (4 " 3e"2) (l - Mcos^ I 


+ 21 


cos'* I ") sin g" - 7e"2 sin^ I" (l-9cos^I") sin 3g'j 


63 


256 e 


M n *) \n' ) ^ (l “ cos I" [^(2 + 33e"2) - (2 - 17e"^) cosM" + 15e"2 sin^ I' 


^K)( 77 ]b^J 


„'2 ^;„2 

2 2 fi _ „» 2 l 2 


a"(l- e"2)^ 


— [^ 2(2 +3e"^) + 15e"^ cos 2g"J 


117 / n ' \ ,42 o ' ^ s in^ I " cos I " 1 

8 (^ng*j ^22 ^..2)7/2 


8 a" ^ (2 + 3e" (l - 3 cos^ I ") + 15e" ^ sin^ I “ cos 


( 181 ) 


cos 2g 


1 .'], 
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I 


dJl_ ^ 

de" de" 


8e 


n \ * 

~J n' ^ a"^ e" j^- (l " 3 cos^ I") +5 sin^ I" cos 2g"J 




fj-e 


" (l - 3 cos2 I") 3 




»2'\S/2 8 J 3 


yu (1 + 4e" sin I " (1 - 5 cos^ I ") 


sin g 


15 . ue 

128 ^ J 


15 , 

256 t’ J 5 


_ e"2)9/2 ( 3-30 cos2 I" +35 cos“ I") 

- 2 (2 + 5e" 2jsin^ I " ^1 - 7 cos^ I "j cos 2g"J 

^(1 - e"" 2 )ii -'~2 [ 2(4 +41e"2 + 18e'"') (l - 14 cos2 I" + 21 cos'* I") sin g" 

~7e"2 (l+2e"2j sin^I"(l-9cos2l"j sin3g"'| 

1^64 -99e"2 + 3(l2-]7e"2j cos2l" 

-*15(2-3e"2) sin2 I" cos 2g"J 


9 /^a 1/ I ^ ^ ^ e"_cos _I" 

128e2 l^n^jVn'J 2 ) 1 /2 ' " 


^ ^2j n ' 2 e" sin2l"cosI" 



22 


(l-e"2)3 


j2 (7 + 3e" + 15 (1 + e " cos 2g"j 


n' 2 e” sin^ I" cos I" 3 , , . 

'2 _ g- 2 j 9/2 ” ¥ '’e a e ^1 “ 3 cos‘‘ I - 5 sin^ i' cos 2g" ) 


4e j 2 


(1 “ 3 cos2 ;i"2 1^2 + (l - e"2 2 ji,^2 j» 


2g" , ( 182 ) 


<91 " dl" 


3 ("El ,2 »2 

^ n ^ a ^ sin 


in I" cos I" 1^- (2 + 3e"2j + 5e"2 cos 2g"J 


_3 , j-i sin I " cos I 

“ 2 *2 J 2 


‘(l-e"2) 


" 2\3/2 8 * 2 ^ J 3 


/<ie" cos I " (11 ~ 15 cos 2 I ") 

a"‘*(l-e"2)s/2 


sin g 


15 , M sin I " cos I " T/ -,\ / ^ » , -I 

32 J 4 gVr^_^„2^7/2 [(2 + 3e 2 ) (3 -7 cos2 I") - 2e“2 (4 - 7 cos2 I") cos 2g"J 
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256 J 5 g"6 /i _ „"2'\9/2 


( 1 -e-)^ 


[2 (4 + 3e"^] (29 - 126 cos^ I " + 105 cos'* I ") sing" 


7e"^sin^I" (7-15cos^I") sin3g"J 


^ \ ^ (l-e"^)*^^sinl" |2 + 33e" ^-3(2-17e"^) cos^I" 

+ 15e" ^ (1 " 3 cos ^ I ") cos 2g"j 


128e^ \n„ / \n 


-L/M/M 2 n'^ sini" (1 -3cos^I") 

32€ P J 22 i, _^«2\2 


( 1 -e'-) 


1 ^ 2(2 + 3e"^) + 15e"^ cos 2g"j 


9 fn' 


4 I ■ * 1 J2^2 


£ / 


1 ' ^ sin I " (1 - 3 cos^ ^ ') 3 r \ 1 

3-2 _p« 2 j 7/2 8 3 “^ sin I" cos I " (2 t 3e"2) + Se"^ cos 2g'J 


_3_ 

4f 


j 2 n^? (1 " 3 cos^ I^) sin I " cos I " -^1 - /3"^ |l * 2 (l - e" ^ cos 2g"^ , (183) 


and 


til , 

■^g" ■ 3 g'' 


1 ^ / 

8 e 


2 o "2 ^"2 ,:„2 T « 


n ^a ‘e ‘sin 


I" sin 2g" - 8 J 3 


/le 


" sin I " (1 - 5 cos ^ I ") 


(1 -e" 2 ) 5 -' 


„ 2 i 5/2 ~ g 


1C /xP “ ^ Sin ^ I " (1 - 7 cos ^ I ") 

A5 , 4 T _ ' 

32 ^ J 4 ,.'5 


_ g» 2^7 /2 S*" 2 g 


15 


fje" sin I " 


256 *^^Js a"6 


2 (4 + 3e " ^) (1 ~ 14 cos^ I" + 21 cos* I ") cos g ' 


- 7e" ^ sin ^ I " ^1 - 9 cos ^ I ") cos 3g"j 


("7 1 n' ^ a"^ e"^ (l - e"^) sin^ I" cos I" sin 2g" 
64e* \n„ y \n / 
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135 

16 ^ Vn*An'j 


n' 2 e"2 sin^ I" cos I" . „ 15 , . „ „ 

J22 ^ ^_g-2)2 sin 2 g - -g- a 2 e 2 5,^2 j- 3,^ 2g 

+ ^ j 2 "c (1 “ 3 cos^ I^)/S"2 j^l + 2 (1 - e"2 )i/2j sir,2 I" sin 2g" 


SUMMARY OF THE DEVELOPMENT 

The short-period terms are given by 


I = r - 


g - g 


r + M 


h = h' 


,9L' 

- A 

-^ = g - Ag 

h' + Ah 


L = 


aH' 

L' + "^7“ = L + AL 


G = G' + 


<?g 


G' + AG 


and 


H ^ H' + 


^2 

(9h 


H' + AH 


Long-period terms are obtained from 


I' = I" 


g ^ g 


h ' - h " 



as; 

as; 

dV 

a.L" 

aL" 

dS^' 

as; 

as; 

dG" 

aG' 

aG' 

as; 

as; 


an ' 

aH 

aH 


= I" + M- 


S + Ag' 


h'- + Ah' 


L' = L" 


( 184 ) 


(185) 
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ds^' aSj' as' 

G' = G" + + ^ + ^ = G" + AG' 

ag' ag' ag' 


and 


as/ as/ as/ 

H' = H" + + + = H" + AH' 

ah' ah' ah' 


(186) 


Secular perturbations and perturbations depending only on g" result from the integration of 
Equations 180: 


and 


M' 


\L - ^ 


J \ aL" 

aL" / 


dt = I" - I, 


Ag" 


• - e 


aF/' 

= g" - g/ 


Ah" = 


!(' 


^ aF/' aF,"' 

n- + + Idt 

aH" aH" 


= h" - h„ 


AL" = 0 , 


AG" 


r/^ ^ ^ 
J\ag" ag" 


dt = G" - G/ 


AH" = 0 . 


(187) 


Then, the perturbations in the Keplerian elements a, e, and I are given by 


a 


a ' 


+ 2a' 


L' 


and 


a ' 



(188) 
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e 


e * + 


1 - e'^ 
e' \L' 


G' / 


_ 1 ~ AG' _ 1 - e"^ MG'\ : 

e" G" 2e"3 \ G" / 

1- ^o"" AG" 


Go" 


and 


e 


0 



(189) 


I = I' + cot I' 


AG AH 
G' H' 


I ' = I " + cot I" 


„ cos I" /AG'\ /AH'\ _ 1 cos I" (l + s in ^ I") ^AG'^ _ 1 ^^^3 

^ V G" H" / ^ sin^ I" \ G" / V H" / ^ sin^ I" \ G' / ^ \H ' 


AG" 

I" = i; + cotig — 


and 


I 


0 



(190) 


POSITION ANO VELOCITY: e ^ 0; I 0°, 180° 

From the elements L, G, H, I, g, h, one can obtain the coordinates and the components of 
velocity as follows. 

Obtain a, e, and I from the following equations: 

L2 

^ ^ M ’ 
e = y 1 - gVL^ , 

and 

I = arc cos ^ (0° < I < 180°) . 
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Now solve Kepler's equation 


E - e sin E - Z 


to obtain E. 

Then compute r from 


Next obtain f from 


r - a(l - e cos E) • 


cos f 


a 

— (cos E - e) 


and 


. ^ _ a G 
sin f - YU ^ • 


Now compute 


and 


Then, 


Aj, - a(cos g cos h - sin g sin h cos I ) , 

G 

- a -jj ( - sin g cos h - cos g sin h cos I ) , 

Ay = a( cos g sin h + sin g cos h cos I ) , 

G 

By - 3 -jj (- sin g sin h + cos g cos h cos I) , 

A = asingsinl , 

Z ® 


B^ - a 17 cos g sin I . 
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Since the system is rotating, 


* A 

rifi. A 

C y 


n* B 

(£ y 




and 


Therefore, 


- n*B 

0 X 




0 
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